arXiv:hep-th/9706222vl 30Junl997 


Confinement in Covariant Gauges 

Martin Schadenf] and Alexander Rozenberg^ 

Physics Department, New York University, 

4 Washington Place, New York, N.Y. 10003 


ABSTRACT 


We examine the weak coupling limit of Euclidean SU (n) gauge theory in covariant 
gauges. Following an earlier suggestion, an equivariant BRST-construction is used to 
dehne the continuum theory on a hnite torus. The equivariant gauge hxing introduces 
constant ghost helds as moduli of the model. We study the parameter- and moduli- space 
perturbatively. For Uf < n quark flavors, the moduli flow to a non-trivial hxed point 
in certain critical covariant gauges and the one-loop effective potential indicates that the 
global SU{n) color symmetry of the gauge hxed model is spontaneously broken to 
Ward identities and renormalization group arguments imply that the longitudinal gauge 
boson propagator at long range is dominated by n(n—1) Goldstone bosons in these critical 
covariant gauges. In the large n limit, we derive a nonlinear integral equation for the 
expectation value of large Wilson loops assuming that the exchange of Goldstone bosons 
dominates the interaction at long range in critical covariant gauges. We hnd numerically 
that the expectation value of large circular Wilson loops decreases exponentially with 
the enclosed area in the absence of dynamical fermions. The gauge invariance of this 
mechanism for conhnement in critical covariant gauges is discussed. 
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1 Introduction 


Some exact results for the measure on the moduli space of supersymmetric theories have 
recently been obtained. Holomorphicity and the exactness of the one-loop /3-function 
greatly restrict the low energy effective action of these models, and one can reach exact 
conclusions about their phase structure. When the number of matter helds is not too 
large, these non-abelian theories conhne]^. 

The absence of fundamental scalars is one of the obstacles which an investigation of 
nonperturbative effects in unbroken ordinary gauge theories has to confront. The moduli 
parameters in this case represent expectation values of composite operators and there 
is apparently no analog to the powerful constraint of holomorphicity for the effective 
potential. It was only recently observed that a translationally invariant quantization of 
ordinary gauge theories on compact manifolds does introduce moduli in a natural way. To 
eliminate normalizable zero modes of the ghosts that would lead to a vanishing partition 
function, the global SU{n) symmetry of the gauge theory is treated equivariantly in a 
quantization which manifestly preserves the isometries of a compact space-time manifold. 
This equivariant BRST quantization [Q introduces a number of global moduli, one of 
which is a scalar in the adjoint representation with vanishing ghost number. A consistent 
covariant gauge hxing procedure thus introduces global moduli parameters even if the 
gauge symmetry is not broken and fundamental scalars are absent at the classical level. 

We here investigate in detail the moduli-space of the covariantly quantized SU{n) 
gauge theory on a torus and study the phase space of its unphysical degrees of freedom 
as a function of the gauge hxing parameters. We consider the most general renormal- 
izable covariant gauges. They depend on three moduli and two gauge parameters: in 
addition to the usual covariant gauge parameter a related to the longitudinal gluon prop¬ 
agator, we also use the freedom to introduce a quartic ghost coupling in the gauge hxing 
functional]^. Although we cannot compute the exact ehective potential, the perturba¬ 
tive analysis reveals hxed points of the parameter and moduli space. These indicate that 
there are at least two asymptotically free phases of the unphysical degrees of freedom for 
Uf < n quark havors. Although no local physical order parameter distinguishes between 
the phases, we argue that a perturbative analysis of the model is consistent only near the 
phase transition. Since the moduli are introduced to avoid infrared problems in the gauge 
hxed action, the measure on this moduli space is gauge dependent. We however argue 
that it rehects non-perturbative properties of the gauge hxed model at hxed points in the 
parameter space. 

The one-loop ehective potential on the moduli space indicates that the phase transition 
is characterized by a condensation of ghost-antighost pairs in the adjoint representation. 
In the inhnite volume limit this leads to a spontaneous breaking of the global SU{n)- 
symmetry to of the equivariantly gauge hxed model. Although this pattern of 

spontaneous symmetry breaking is reminiscent of monopole condensation, the composite 
local scalar held in this case is gauge dependent and not a physical order parameter. One 
has to consider non-local Wilson |^- or equivalently, magnetic]^]- loop operators as gauge 
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invariant order parameters for confinement. 

The scale of symmetry breaking is found to be RG-invariant and proportional to A asp 
on a line in the parameter space. For a fixed number Uf < n fermions in the fundamental 
representation, there are two hxed points on this line which merge for Uf = n. The 
corresponding covariant gauges we call critical covariant gauges (CCG), since perturbation 
theory appears to be selfconsistent in these gauges also for unphysical degrees of freedom. 

Ward identities are used to show that the correlator of the composite local operator 
oc with the conserved color current in CGG has a Goldstone pole in 

the thermodynamic limit. We infer from this, that {d-A{x) (9-R(0)) also has Goldstone 
behavior at large distances. This behavior is shown to be compatible with the Ward 
identities in CCG. In these gauges with quartic ghost interactions. Ward identities do not 
completely determine the longitudinal gluon propagator, but rather relate it to correlation 
functions in the ghost-antighost channel. 

The contribution of the long-range longitudinal interaction to the expectation value 
of large Wilson loops is then obtained in leading order of n. Numerical solution of the 
non-linear integro-differential equation we derive in this approximation shows that the 
expectation value of large circular Wilson loops decreases exponentially with the enclosed 
area over 12 orders of magnitude. 

The structure of this article is as follows. In section 2 we obtain the continuum action 
of the equivariantly quantized theory in general covariant gauges on a torus. This step 
is crucial, since it gauge fixes also the long wavelength excitations of the Yang-Mills 
theory without destroying the homogeneity and isotropy of the model. The equivariant 
construction introduces global fields, which we subsequently treat as moduli. In section 3 
we study the renormalization flow behavior of the gauge parameters in such generalized 
covariant gauges. We locate the ultra-violet fixed point of the gauge parameter a, which is 
related to the longitudinal gluon propagator at large momentum transfers. In section 4 we 
determine the non-trivial hxed points of the moduli-space introduced by the equivariant 
gauge hxing. We hud that the minimum of the effective potential is at a non-trivial value 
for one of these moduli and that the scale of this minimum is proportional to A asp in 
CGG. The critical exponent of the non-trivial moduli vanishes in these critical gauges. 
The one-loop effective potential also indicates that the global SU (n)-symmetry is broken 
to in critical gauges. In section 5 we exploit the consequences of this symmetry 

breakdown using Ward identities of the model. These provide the connection between the 
perturbative indication that the SU (n)-symmetry of the model is spontaneously broken 
and long-range effects. We establish the existence of n{n — 1) unphysical Goldstone 
bosons associated with the global symmetry breaking by studying the (^“-correlator with 
the conserved color current. Due to mixing between (9“ and d-Ahy renormalization, 
the Goldstone bosons also dominate the long-range behavior of the longitudinal gluon- 
propagator in CCG. We show that this behavior is consistent with the Ward identities. 
In section 6 we use these results to approximately evaluate the expectation value of large 
circular Wilson loops in leading order of n. We derive a non-linear integral equation which 
resums the longitudinal interaction in leading order of n and we solve it numerically. 
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Section 7 is a summary and discussion of the results. 


2 SU{n) gauge theory in covariant gauges 


To investigate phases of a quantum field theory, one has to control it in the thermodynamic 
limit of inhnite space-time volume. We will study order parameters of an Euclidean SU{n) 
gauge theory on a torus as the space-time volume of the torus becomes large. 

A lattice regularization might appear to be the only reasonable dehnition of a gauge 
theory, since it does not depend on perturbation theory. Although this is a perfectly 
regularized statistical theory on a point manifold, the physical volume of any hnite lattice 
shrinks exponentially as the coupling constant is tuned to the critical value. There is 
little hope of studying long-range effects numerically in the critical regime, and current 
simulations are restricted to a few fm^ of physical volume. These numerical investigations 
suggest conhnement, but deconhnement at a couple of GeV excitation energy cannot be 
excluded by these studies. Experimental evidence for the absence of colored (and/or frac¬ 
tionally charged) asymptotic states is considerably more restrictive and likely to remain 
so for some time. 

A theoretical investigation of the thermodynamic limit of lattice gauge theory near its 
critical point is rather complicated. It has been pointed out|p, that the Gribov problem[0] 
cannot be ignored in a study of the lattice in this critical limit. A proposed effective 
description of lattice gauge theory at weak coupling introduces an additional parameter P, 
which effectively constrains the conhguration space to the fundamental modular region in 
the thermodynamic limit. It was shownthat a non-vanishing value for this parameter 
spontaneously breaks the BRST-symmetry of a lattice in Landau gauge^. 

We will not attempt to derive the thermodynamic limit of a lattice gauge theory here, 
but propose to investigate instead the consistency of perturbation theory in a covariantly 
gauge hxed continuum theory with presumably the same critical behavior as a lattice 
gauge theory. Let us hrst perturbatively dehne the continuum gauge theory on a compact 
space-time manifold such as a hypertorus. In order to effectively exploit the space-time 
symmetries of the thermodynamic limit, we consider only gauges which manifestly pre¬ 
serve all the isometries of the torus, i.e. general covariant gauges. The torus with periodic 
boundary conditions for the gauge and ghost fields is chosen for the following reasons: 

• analytic continuation in the dimension D of the hypertorus To of perturbative 
expressions is straightforward. The perturbative dimensional regularization does 
not break the gauge- (or better BRST-) invariance explicitly. It allows one to 
efficiently extract the ultraviolet finite renormalized quantities in the limit D ^ A 
and their critical exponents. 

^The BRST-symmetry of a covariantly fixed finite lattice has to be broken, since the Euler character 
of the lattice gauge group vanishes and the gauge fixed partition function of any BRST-symmetric hnite 
lattice model would therefore vanish as welllpi. 
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• The space of gauge orbits of an SU (n) theory on a hypertorus with periodic bound¬ 
ary conditions is connected. We thus only have to consider perturbation theory 
around a single classical gauge orbit. The absence of nontrivial topological sectors 
on the finite torus with periodic boundary conditions might oversimplify our study 
of the thermodynamic limit of this model. The non-perturbative effects we wish 
to investigate, however seem to be more easily described in a toroidal geometry of 
Euclidean space-time. 

• An SU{n) gauge theory defined on a torus resolves the topological ambiguities of 
the covariant gauge fixing procedure we employ pj. These topological anomalies are 
present only on base manifolds with non-trivial 3-cycles. 

• there is no topological obstruction to introducing fermions [pT|] in the fundamental 
representation on a torus with periodic boundary conditions. This would not be 
possible with twisted boundary conditions for the gauge field^. 

• The torus is a flat Euclidean space-time manifold. Translations and the associated 
Fourier - analysis are straightforward, and an explicit mode expansion is manage¬ 
able. 

These topological properties of a torus with periodic boundary conditions facilitate the 
perturbative analysis considerably. We will repeatedly take advantage of the translational 
and rotational invariance of gauge-dependent correlation functions in the thermodynamic 
limit of inhnite space-time volume in covariant gauges. 


2.1 Equivariant gauge fixing 

Thermodynamic arguments for gauge dependent correlation functions are valid in co¬ 
variant gauges, which by definition preserve the isometries of the base manifold. The 
Goldstone theorem in this case also gives information about the long-range behavior of 
gauge-dependent Green’s functions. 

It is well known that harmonic ghost modes pose an obstacle to the definition of a 
covariantly gauge-fixed partition function on compact manifolds ||11||. These zero-modes 
of the Faddeev-Popov ghosts are a consequence of the invariance of covariant gauges 
with respect to global SU{n) rotations. The zero modes are conventionally removed by 
fixing the global SU{n) invariance at a particular space-time point |^. Although this 
procedure preserves the BRST-algebra, such pointed gauges single out a space-time point, 
and thus break the translational invariance of gauge-dependent correlators explicitly. If an 
unbroken SU{n) gauge theory indeed confines color, the usual argument that this point- 
defect can be disregarded in the thermodynamic limit of gauge dependent correlation 
functions does not necessarily hold, since color correlations must be strong and of long 
range. It is therefore not obvious that translational invariance of gauge dependent Green’s 
functions, such as gluon- or ghost-propagators, is restored in the thermodynamic limit of 
pointed gauges. If translational and Euclidean rotational invariance of gauge dependent 
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correlators on the other hand cannot be explicitly maintained, the motivation for choosing 
a covariant gauge becomes questionable. 

Recently a solution to the problem with harmonic ghost zero modes was proposed, 
which manifestly preserves the invariance of the gauge hxed theory with respect to all 
isometries of a compact manifold|Q. The global S'17(n)-invariance is treated equivariantly 
in this construction and the usual held content of an SU{n) gauge theory augmented 
by constant global ghosts whose canonical dimensions and ghost numbers are given in 
Table 1. 


held 

A^,{x) 


c{x) 

c(x) 

h{x) 


UJ 

a 

cr 

7 

7 

dim 

1 

3/2 

0 

2 

2 

0 

0 

4 

4 

2 

2 

0n 

0 

0 

1 

-1 

0 

2 

1 

-2 

-1 

0 

1 


Table 1. Dimensions and ghost numbers of the fields. 


The action of the nilpotent BRST 

sA^{x) = 
sc(x) = 

SUJ = 
S(j) = 

sc{x) = 
sb{x) = 
sa = 
sa = 

S7 = 

sy = 
S'lpi^x) = 
si’iix) = 


operator s on the helds of Table 1 isQ 

D ^ c ( x ) - [ u , A ^{ x )] 

-[ u , c { x )] - |[c(a:),c(x)] - -0 
^ 9 

-^[o^,a;] + 0 

-[a;,0] 

-[a;, c { x )] + h { x ) 

-[ uj , h { x )] + [0,c(a;)] 

-[a;, cr] + d 
-[ uj , a ] + [0,cr] 

-[^,7] +7 
-[a;, 7] + [0,7] 

- uji > i { x ) - gc { x ) iii { x ) 

- iji { x)u - g ' ipi { x ) c { x ) 


where 

D^c{x) = d^c{x) + g[A^{x),c{x)] 
is the usual covariant derivative for the adjoint representation. 


( 2 . 1 ) 

( 2 . 2 ) 


It is straightforward to show that the BRST-operator dehned above is nilpotent on 
any element of the graded algebra constructed from the helds of Table 1: 

= 0 (2.3) 

^We will often suppress color indices. Except for 0i, 0i, vectors of the two fundamental representations 
of SU{n), all fields are traceless, anti-hermitian nx n matrices in this notation, and transform under the 
adjoint representation. 
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The constant ghost a; generates global gange transformations of all the fields in the BRST 
algebra, except itself. One therefore can restrict observables to the eqnivariant cohomology 

S, 

S = {C>: 1^ = 0; sC> = 0,0^ (2.4) 

where T is itself (^-independent. Since one is interested only in expectation valnes of 
gange invariant fnnctionals of A, and the notion of physical observables can be fnrther 
sharpened to fnnctionals in the eqnivariant cohomology with vanishing ghost nnmberQ 
The most general power connting renormalizable classical action S with the symmetry 
O depends on two gange parameters a and (3 = a5, 

S = Sc[g, rrii] A, -ipi, -0^] - 1 - sWcrig, a, /?, 7 , a; A, c, c] (2.5) 

where Sc is the Yang-Mills action of an SU{n) gange theory with nj qnark flavors in the 
fnndamental representation 

Sc= f -^:rF^,{A)F^M) + ^^/J0 + m,)^/J, (2.6) 

■Jt Z 

The Enclidean Dirac operator here i^ 

Ipipi = (2.7) 

and the field strength is related to the connection by 


F^v{A) — d^Ay — duA^ -|- g[A^^ A^] 


(2.8) 


Note that a CP-breaking term proportional to the Pontryagin nnmber Tr F AF vanishes 
on a torns with periodic bonndary conditions^ We therefore do not inclnde snch a term 
in (|2.6|) . The discrete CP-symmetry in this sense is not broken on any finite torns. 


The BRST-exact gange fixing term of the action (|2.5|) was obtained in 

^GF = 2 


IT 


dx Tr [d^c{x)A^{x) - ac{x)h{x)—a5gc{x)c{x)c{x) + —^c{x) + —ac{x)] (2.9) 

2 g g 


For a > 0 the Nakanishi-Lantrnp field h{x) in ( p.5|) can be eliminated nsing the eqnation 
of motion 

h{x) = iA/g - d-A{x))/a - 5g[c{x), c(x)], (2.10) 


since observables do not depend on h{x). 

®To avoid global anomalies due to large gauge transformations, one could consider only local observ¬ 
ables. However, the Chern-Simons term appears to be the only gobal physical observable with such a 
gauge- ambiguity]^. The restriction to local observables is therefore not necessary on a torus. 

®The Euclidean Dirac matrices 7 ^ satisfy 7 ^ 7 ,^ -I- = 2(5^,^!. 

^ “Toron” sectors with non-vanishing Pontryagin number exist on a torus with twisted boundary condi¬ 
tions, which however can only be imposed in the absence of quarks in the fundamental representation]^. 
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A covariantly quantized SU (n) gauge theory on a torus is thus described by the tree- 
level action 


*S'o = Sc + “2 dx Tr 


IT 


2a 


{d-A{x)Y -I- 6c{x) D^-dc{x) + (1 — 5)c{x)d-D'^c{x)-\- 


-|-a5(l — 5)g^c{x)c{x)c{x)c{x) — S[c{x), c(a ;)]7 -|- 


^,T-—7d-A- 

2ag^ ag 


-ac{x)c{x) — a{l — S)(j)c{x)c{x) - -acf) -\ — 'jc{x) -\ — ac{x) 

9 9 9 


( 2 . 11 ) 


The constant Grassmann ghosts 7 and a in ( |2.11| ) implement the constraints 

J dx c{x) = 0 

J^dxc{x) = 0 ( 2 . 12 ) 

which remove the (on a finite torus normalizable) constant modes of the ghost-helds in 
a translationally invariant fashion. In the framework of TQFT one can show that the 
partition function of the gauge-hxed quantum theory on a compact space-time manifold 
otherwise would vanish due to these zero-modes of the ghosts ||TI|, ||, |[. The surface term 
■ A in ( p.ll| ) vanishes on a hnite torus with periodic boundary conditions. The 
absence of this term on a torus will eventually give rise to a non-trivial effective potential 
on the moduli-space. We retain this surface term in order to manifestly exhibit the 
BRST-symmetry of ( p.ll|) . 

The action (^. 11 |) is formally invariant under the action of s dehned by: 


sA^{x) = D^c{x) 

sc{x) = {f/g — d-A{x))/a — 6g[c{x),c{x)] 
sa = d 
57 = 7 

S'lpiix) = -gc{x)'ip{x)i 


sc(x) = -llc(x),c(x)] - 
sf = 0 
sd = [0,cr] 

_ 57 = [0,7] 
s-ifiix) = -g'fi(x)c(x) 


(2.13) 

which is on-shell nilpotent only on the set of globally SU(n) invariant functionals. Apart 
from an unorthodox transformation of the anti-ghost c and the presence of additional 
constant ghosts, s is similar to the conventional BRST-transformation, and dehnes the 
equivariant cohomology of physical observables. 

The constraints ( | 2 . 12 |) are invariant with respect to (|2.13| ) at points where 

J dxsc{x) = 0 and J dxsc{x) = 0 (2-14) 

Dividing by the hnite volume of the torus and taking the thermodynamic limit, (|2.14|) 
relates the global bosonic ghosts 7 , 0 to vacuum expectation values of composite dynamical 
helds 

( 7 ) = c^Sg^ ([c(a;), c(x)]) and (0) = -g"^ {c^{x)) (2.15) 
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The relations ( [^.151 ) show that the equivariant BRST-symmetry of the model is unbroken 
in the thermodynamic limit only at very specific values for the global ghosts. In what 
follows, we will see that the first of these relations in (p.l5|) is satisfied non-trivially in 
certain gauges and will explore the consequences. 


The expectation value of an observable (9 G S is formally given by the path integral 

{O) = A/” y d(f)d(jdad^d'^ jJ [V'l/jVipVAVcVd] O expS'o, (2.16) 

which defines the perturbative loop expansion. We show in Appendix A that (|2.16|) is 
generally normalizable, i.e. that ( 1 ) ^ 0 for an SU{2) group. The global bosonic ghosts 
0, a and 7 introduced by the equivariant BRST-algebra are moduli parameters of the 
theory. We indicate in (|2.16|) that the integration over this finite dimensional moduli- 
space would usually be performed after the path integral over the dynamical fields. As 
far as the dynamical fields are concerned, the moduli are parameters of the action. Note 
that the moduli space possesses a trivial BRST-symmetry. 

With finite fermion-masses m*, and/or anti-periodic boundary conditions for the fields 
'ipi and fji, only the constant modes of the connection A^(x) are potentially troublesome 
in a perturbative treatment of the theory. One can eliminate these constant modes by 
requiring twisted boundary conditions p. These boundary conditions however break the 
global SU (n) invariance of the gauge theory explicitly^] and effectively introduce external 
color flux|^. Although this is arguably the most elegant treatment of the purely gluonic 
theory on a torus, twisted boundary conditions require the absence of fields in the fun¬ 
damental representation. They also complicate the mode expansion and thus the infinite 
volume limit, especially since the space of gauge orbits in this case is disconnected. The 
global SU{n) symmetry of the equivariantly gauge fixed theory is, as we shall see, sponta¬ 
neously broken even for periodic boundary conditions at certain points in the parameter 
space. The theory is thus seen to react extremely sensitively to any external color flux. 


The constant part of a field configuration A^(a;) satisfying periodic boundary con¬ 
ditions 


Af,{x) = -a^ -h Af,{x) with 


Jt 


dx Afj_{x) = 0 


(2.17) 


can be treated as another moduli on which the perturbative expansion depends. Since 
we only evaluate correlation functions of local operators, the integration over this moduli 
space is governed by the corresponding effective action. Due to translational invariance of 
the constant modes, this effective action is proportional to the volume of space-time. In 
the infinite volume limit only the immediate vicinity of the absolute minimum (possibly 
degenerate) of the effective potential therefore contributes to the integration over the 
moduli-space. Contrary to supersymmetric gauge theories, we do not obtain the complete 
dependence of the effective potential on the moduli parameters. The assumption that 
perturbation theory is a good asymptotic expansion in conjunction with some physical 
criteria, however greatly constrains the minimum of the effective potential. Only the 


®the equivariant construction is then in fact not necessary, since twisted boundary conditions also 
eliminate the constant ghost modes 









minima which are stable fixed points can be consistently incorporated in a perturbative 
expansion of the theory, since moduli tend to flow to such fixed points as the coupling g 
becomes critical. 


From a physical point of view, the minimum of the effective potential should occur for 
= 0 , since translational invariance implies that 

{Af,{x)) = ^ {ttf,). (2.18) 

The discrete rotation group of a symmetric torus would be broken unless all (a^), fi = 
1,..., H are equal. Equivariant observables such as Tr F^yFp„ or furthermore 

would register a breakdown of the S'(9(4)-symmetry in the continuum limit for (Ap) ^ 0, 
unless this phase is gauge-equivalent to one with vanishing = 0 . A phase with 

nontrivial (A^) is therefore acceptable only if 

(a^) = U\x)dpU[x) (2.19) 


for some U{x) G SU{n) on the torus. In this case the theory is however physically 
equivalent to the one with (a^) = 0. Since the torus has only 1-cycles, (|2.19|) implies that 
{ttp) is abelian, i.e. [(a^), {a^)] = 0. In Appendix A we find these abelian Gribov copies 
of the trivial vacuum Ap = 0. They are not removed by the equivariant construction on a 
torus. The contribution to {Ap) from these gauge copies however cancel. Gauge invariant 
local physical observables do not depend on which of these copies is selected, and we may 
as well perturbatively expand around the trivial configuration (Ap) = 0. 


The global SU{n) symmetry is also spontaneously broken by a non-vanishing expecta¬ 
tion value of any of the constant bosonic ghosts (/>, a or 7 of the equivariant gauge fixing. 
None of the local functionals in the equivariant cohomology S are order parameters for 
a spontaneous breakdown of the global SU{n) symmetry. We will see, the the perturba¬ 
tive evaluation of correlation functions in certain gauges is however sensitive to whether 
and how the global S't/(n)-symmetry is broken. The situation is analogous to ordinary 
spontaneous symmetry breaking if one restricts the observables of the model to group 
invariants: the perturbative calculation of “invariant” correlations depends on whether 
the vacuum breaks the symmetry or not, but none of these invariant observables indi¬ 
cates which of the degenerate symmetry breaking vacua has been selected. Furthermore, 
the Goldstone bosons of the symmetry breaking do not lead to massless singularities in 
correlations of the invariant “observables”. For such a restricted set of observables, the 
symmetry breaking is in a sense hidden and its effects appear to have a dynamical origin. 


Non-vanishing (0) or {a) implies that the ghost number is not conserved. The expec¬ 
tation value of the bosonic ghost 7 however conserves the ghost charge and there is no a 
priori reason for it to vanish in the equivariant theory. 
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3 Fixed points of the gauge parameters 


We begin by investigating the critical behavior of the dimensionless gange parameters a 
and 6 for g ^ 0. For 5 = 0, the gange parameter a was observed to flow to a non-trivial 
hxed point 1^. We will see that a flows to a hxed point a^o for a wide range of valnes of 
the gange parameter 6. 

Mnltiplicative renormalization constants relating bare (B) to renormalized (R) dy¬ 
namical helds and parameters are dehned as usnal 

V’h) = = Ab = {cBiCb) ^2 

ZqQr = Qb Z^oiB = o^B ZsSr = 5b 


At the one loop level, 5 does not renormalize [^. Connter terms for the glnon polarization 
determine as well as Z^. Note that the usnal Slavnov-Taylor identity for the longitu¬ 
dinal part of the gluon propagator]^ only requires equality of Z^ and Z 3 for 5 = 0. In 
generic covariant gauges the equation of motion of the Nakanishi-Lautrup held ( p.lOl ) also 
involves the composite operator [c(x),c(a:)] and Z^ 7 ^ Z^. To one loop this is conhrmed 
by an explicit calculation of the renormalization constants of (p.l|) . In the MS-scheme 
one obtains in = 4 — 2^ space-time dimensions 


^2 

^3 

^3 

Zql 

Z, 

Zs 


1 - 

IdiT^ne 

- + 0{g^) 




1 + 


”(3 

«) + 0{g^) 




1 + 

167r^£ 

(f" 

- 1"/ - 

) + 0(g 



1 + 


(f" 

- 1^/ - W 

+ 5(1- 

5)Q;n 

1 + 0 ( 9 “) 

1 - 

32TT^e 

l3, + Oig^) = l- 

( 

— 

D 

in/) + 0{g^ 

1 + 

0{9^) 







(3.2) 


(The fact that Z^ = Z^ also for 5 = 1 is not a consequence of the usual Slavnov-Taylor 
identity [0].) Note that the renormalization constant Z 3 of the dynamical ghosts does 
not depend on 5 at the one-loop level. In ( ^.21 ) and most of the following we dropped the 
index (R) for renormalized quantities. 


For a = Q!oo 


(^) 


13n — 4n/ 
3n(l-25 + 252) 


(3.3) 


Zq is independent of the renormalization scale fi to order g"^. It is readily seen that a^o is 
the stable hxed point of the gauge parameter a which determines the asymptotic behavior 
of gauge dependent correlations. 
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The relation (|3.1|) between bare and renormalized gauge parameter implies for the 
renormalized a 

+ = 0 (3.4) 

d In p d In /i 


in any gauge other than Landau gauges. Integration of (|3.4|) using the 1-loop expression 
for Zn and 


dg 


= Pig) = 


g 


11 2 

-Po + ... with I3q = -n - -Uf 


(3.5) 


dln/i Ibvr^ 

shows that the gauge parameter a approaches the hxed point a^o of ( p.3|) for sufficiently 
small coupling g as 

13ti—4ti ^ 

a — a^o = (1 + ) (3.6) 

where C is an integration constant related to the choice of gauge at a particular (small) 
value of the coupling. C is related to a in a similar fashion as Is .asp is related to the 
coupling g: it is the RG-invariant “gauge-parameter” which governs the asymptotic flow 
of a as /i —cx). Note that the critical exponent in (|3.6|) is independent of S and positive 
in the asymptotically free regime when 
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(3.7) 


This implies that Landau gauge is not the asymptotically stable hxed point in the gauge 
parameter space if the number of families is not too large. Landau gauge becomes the 
stable hxed point only when the exponent in (|3.6|) is negative, i.e. for > nj > 


havors. In the case of QCD this occurs only for 5 — 8 families. Note also that a tends to 
ttoo slower than The retained term of Za therefore dominates over the neglected 

terms of order g^ in (|3.4|) at sufficiently small coupling. For Uf = n, a approaches ctoo 


proportional to g. We remark that the hxed point a^o as well as the critical exponent 
in ( |3.6|) tend to hnite, non-vanishing values in the limit of large n, a case which will 
interest us. In general, describes a surface in the space of continuous parameters 

a,6, g,nf/n of the model. We include the ratio nj/n as a parameter of the model, since 
asymptotic freedom for instance depends on its value. 

For hxed Uf/n, the perturbative expansion of an SU{n) gauge theory is selfconsistent 
in the vicinity of the curve in the coupling space g, a, 6 corresponding to (7 = 0 , and 
a = acx)(^) given by ( p.3|) . The coupling constants g and a can be traded for the RG- 
invariants Aasp and C. RG-invariant quantities only depend on Aasp, C and 6, and 
physical gauge invariant correlations should furthermore not depend on the last two. An 
immediate consequence of this analysis is that the leading deep Euclidean behavior of 
the longitudinal gluon propagator is ehectively described by the tree-level propagator 
in the gauge aoo{S): corrections to this asymptotic behavior vanish logarithmically for 
sufficiently large momentum transfers. At the hxed point aoo the term d/da in the RG- 
equation only induces corrections that are analytic in g and which vanish for 77 ^ 0 . 

It is perhaps surprising that a gauge parameter such as a depends on the renormal¬ 
ization scale at all. a was introduced as a coupling constant of the TQFT on the gauge 
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group [Q and the scaling-functions of a TQFT should vanish. The coupling constants of 
a TQFT are therefore not expected to depend on a renormalization scale. We do not 
have an exhaustive explanation for this apparent paradox. However, if one distinguishes 
between the dehnition of the TQFT and the perturbative calculation giving the apparent 
scale dependence of a, this may not be so paradoxical. We wish to point out that only 
one Gribov copy of a gauge held conhguration contributes in the perturbative evaluation 
of a Green’s function (the one in the vicinity of A^{x) = 0 ). On the other hand, the 
conclusion that a parameter of the TQFT does not depend continuously on a scale re¬ 
quires the contributions from all stationary points, that is from all Gribov copies of a 
conhguration. It is therefore conceivable that the asymptotic behavior of the longitudinal 
gluon propagator evaluated within a single Gribov region depends on the renormalization 
scale, whereas it would be scale independent, if the contribution from all other Gribov 
copies were considered as well (it is known that the “size” of a particular Gribov region 
is scale dependent 0). Perturbatively one verihes only that gauge-dependent quantities 
evaluated in a single Gribov region, asymptotically require an adjustment of the gauge 
parameter with the scale /x. At the hxed point Ooc this scale dependence vanishes and 
the contribution from other Gribov copies to the asymptotic behavior of gauge dependent 
Green’s functions can be ignored. We hnd the existence of such a hnite hxed point ctoo 
remarkable. 

The gauge invariant correlation functions one is ultimately interested in should not 
depend on whether all gauge copies of a conhguration are taken into account or only 
one. Since the “number” of such copies does not change within a connected sector of the 
orbit space 1^, the diherence between the two evaluations of gauge independent Green’s 
functions within a single topological sector can be absorbed in the normalization of the 
path-integral. The orbit space on a torus is connected, and one should therefore obtain 
the asymptotic expansion of physical correlation functions by perturbing around a single 
classical solution satisfying the gauge condition. (Seep for a discussion of the dependence 
of the normalization on the sector for an SU{2) gauge theory on S 4 .) 

We have found the hxed points in the parameter space. In the next section we consider 
hxed points of the moduli space. We will see that perturbation theory is consistent in a 
phase where the global SU{n) symmetry is spontaneously broken to for particular 
values of the gauge parameter 6 . Although this consistency requirement for the perturba¬ 
tive expansion is gauge dependent, the fact that the global SU{n) symmetry is broken to 
perhaps is not. We will hnd that our moduli space is however consistent with such 
a symmetry breaking only on a subset of the hxed points in the parameter space. For 
other values of the parameters, non-perturbative conhgurations are probably relevant even 
asymptotically - just as contributions from other Gribov regions probably would stabilize 
a in gauges a 7 ^ ttoo- At particular points in the parameter space, perturbation theory 
indicates that the moduli space might be adequate to describe this symmetry breakdown 
and that other non-perturbative phenomena could be of sub-leading importance. 
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4 Fixed points of the moduli space 


Perturbation theory is not entirely determined by the UV-£xed points of the parameters. 
It also depends on the phase, which usually is encoded in the expectation values of scalar 
helds. In the case of gauge theories these are composite operators, whose expectation 
values generally will not be described by the moduli of our model. We can however per- 
turbatively examine whether this restricted moduli space allows for non-trivial expectation 
values at particular values of the model parameters. 

Perturbation theory on the torus in covariant gauges depends on the moduli-space M. 

= {y, 0, a, a^} (4.1) 

The measure on this moduli-space is induced by integrating out the dynamical helds, 
which gives the effective potential for the moduli. In TQFT’s the semi-classical measure 
on the moduli-space is exact, and in supersymmetric theories it is also severely constrained 
by holomorphicity|l[]. 

In contrast to supersymmetric gauge theories, this measure on the moduli space will 
not be obtained exactly for the gauge theories that interest us. In a translationally in¬ 
variant model, the effective action on the moduli space is however proportional to the 
space-time volume of the manifold. In the thermodynamic limit the integration over the 
moduli-space is thus effectively restricted to regions arbitrary close to absolute minima of 
the effective potential. To evaluate the integral over the moduli space in the thermody¬ 
namic limit in four space-time dimensions, we in principle only need to know the scale of 
these minima and the behavior of the effective potential in their immediate vicinity. If a 
compact group relates degenerate absolute minima of the potential and the observables 
are invariants under this group, it furthermore suffices to evaluate the moduli at one of 
the degenerate absolute minima only. The group is then spontaneously broken to the 
subgroup which leaves this minimum of the effective potential invariant. 


Suppose for a moment that the effective potential of the asymptotically free theory in 
four dimensions has a unique RG-invariant non-trivial absolute minimum modulo group 
invariance. In the absence of dimensionful parameters, the scale k associated with this 
minimum|^ is a function of the renormalization point /i and the dimensionless coupling 
constants g, a, 6 only. Since the minimum is a hxed point of the moduli space and the 
theory is asymptotically free, n must be proportional to Kasp and its dependence on /i 
and g at weak coupling therefore is 


In 


K 




Ibvr^ 


A, 

A Ibvr^ 


+ Ck -|- 0{g‘^) 


where 


A — 
A = 


11 


n — 




¥f 

- fnnf + 


III 

n 


( 4 . 2 ) 


( 4 . 3 ) 


®The canonical dimension of (j) in (O) vanishes, but the effective potential conserves ghost number 
(although its minima need not to), and is thus a function of acj) rather than (j) itself. 
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are the first two coefficients of the /3-function. Since the leading coefficient oil/ in (^4.2|) 
cannot vanish, k has to be the scale of the one-loop effective potential also. In fact, k has 
to be associated with a non-trivial absolute minimum of the one-loop potential in the case 
we are considering, since the RG-invariant absolute minimum of the effective potential 
defines a mass scale. Conversely, a non-trivial absolute minimum of the one-loop effective 
potential which is RG-invariant and unique modulo group invariance defines an asymptotic 


mass scale which can be related to other definitions of A^sp by a two-loop calculation|15 


It should be emphasized that an RG-invariant minimum of the effective potential on the 
moduli space gives an asymptotic scale, whose coefficients in the expansion (|4.2|) can be 
reliably calculated perturbatively. 


The coefficient of the term in (f4.2|) is determined by the one-loop effective potential 
and must be —(47r)^//3o at a fixed point of the moduli space. Depending on whether 
A ASP n, or A asp -C k at weak coupling, one otherwise concludes that the only fixed 
point of the moduli space is the trivial one with k/Aasp = 0 or the one at infinity. In 
the latter case the perturbative evaluation of the potential is certainly inconsistent, since 
it gives rise to an arbitrary large mass scalej^. It should be emphasized that this first 
coefficient of the expansion ( [4.2|) is entirely determined by the classical term of the effective 
potential and its divergent one-loop contribution. Since the effective potential is only UV- 
divergent, a perturbative calculation of this coefficient is reliable in asymptotically free 
gauge theories. 


Whether k is truly marginal for > 0 is determined by the coefficient of the In g"^ term 
in (|4.2|) . It is given by the one-loop anomalous dimension of the order parameter. Since 
the effective potential is itself RG-invariant, the coefficients of the 1/g"^- and In^f^-terms 
in fact are related by the RG-equation: near a fixed point in the parameter space, k is of 
the form ( [4.2| ) only if the critical exponent of the moduli parameter vanishes. The relation 
between the first two coefficients in (^) also greatly facilitates the search for fixed points 
of the whole moduli space, since one merely needs to compute the critical exponents of 
the moduli on the curve (|3.3|) in the gauge parameter space. As shown in Appendix B, 
the anomalous dimension of a moduli is also determined by the divergence of the one-loop 
effective potential and the associated classical counterterm. 


A two-loop calculation of the effective potential would in addition uniquely determine 
the coefficient in ([4.2|) and thus relate k to any other asymptotic scale parameter 
quantitatively!^. To qualitatively analyze the moduli space, the relation of k to other 


asymptotic scale parameters is not required and will not be obtained here. We neverthe¬ 
less emphasize, that a unique RG-invariant non-trivial absolute minimum of the effective 
potential on the moduli space is perturbatively related to other definitions of the asymp¬ 
totic scale parameter in ordinary gauge theories. 


It might appear that the converse is not necessarily true, and that a non-trivial abso¬ 
lute minimum of the one- and two-loop effective potential with the above properties does 
not necessarily imply the same for the full effective potential. The argument above how- 

^°the usual A^sp should set the mass scale and the asymptotic expansion is not expected to depend 
on masses much larger than A^sp- 
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ever shows that this would question the validity of perturbative results in ordinary gauge 
theories since non-perturbative effects would have to redehne the perturbative asymptotic 
sea leQ. Perturbation theory by itself ought to be adequate in determining critical expo¬ 
nents etc. of ordinary gauge theories, and the perturbative asymptotic scale can be given 
a physical interpretation in this case. This implies that the RG-invariant non-trivial abso¬ 
lute minimum of the two-loop effective potential which dehnes an asymptotic scale cannot 
simply be absent from the full effective potential. This argument for the existence of such 
a scale-dehning absolute minimum also in the full effective potential does not mean that 
the two-loop effective potential is exact in ordinary gauge theories. It only implies that a 
RG-invariant scale set by the absolute minimum of the effective potential can be used to 
normalize the model also perturbatively and that the full effective potential must have a 
non-trivial absolute minimum, if the one-loop potential does and if this minimum is RG 
invariant. In the following we therefore only verify whether the one-loop effective poten¬ 
tial on the moduli space has non-trivial RG-invariant absolute minima at hxed points of 
the parameter space. 


Let us hrst assume that ghost number and angular momentum are conserved, i.e. that 
(0) = {a) = (a^) = 0. In this case we need to consider only the dependence of the effective 
potential on the moduli 7 . Generalizing the calculation for an SU{2) gauge group in[@] 
to SU{n ), we obtain in Appendix B the dependence of the one-loop effective potential on 
the eigenvalues Xi,i = 1 ,..., n of the traceless n x n hermitian matrix iy 


Gl— ioop ^ ' 

l<2<j<n 


(A, - 

327r2 


In 


(A, - X,f6^ 




(4.4) 


Note that the full effective potential also depends only on eigenvalue differences, and that 
the logarithm in ( |4.4|) is typical for the one-loop approximation. Grucial for the following 
is that (^4.4|) has a unique and non-trivial absolute minimum (modulo SU{n)) and the 
dependence of the associated scale k in (|4.4|) on the renormalization point pi and the 
coupling constants. To this order in the loop expansion we obtain in Appendix B that. 


Anpi^ 


(4vr)^ 

n6‘^ag‘^ 


+ 1 -7e + 0(hi/,/) 


( 4 . 5 ) 


in the MS-scheme. Gomparing with ( |4.2|) we see that k is an asymptotic scale only in 
gauges where 

^ /^o (4.6) 

The intersections of this curve with aoo((5) of ( p.3|) occur for gauge parameters 5 = 5^ 


lln — 2nf i: \j2{n — nf)(lln — 2n/) 


D — 2 space-time dimensions the order parameter is dimensionless and the effective potential does 
not define an asymptotic scale and there is no spontaneous symmetry breakingp^. In supersymmetric 
gauge theories in four dimensions, non-perturbative corrections are asymptotically important [T^ and 
their contribution is required to unambiguously define an asymptotic scale. In theories which are not 
asymptotically free, such as 4>^, our argument obviously does not apply. 


( 4 . 7 ) 


15 












For a given value of Uf/n, the two points in the gauge parameter space (Q;oo(^i,), and 
(aoo(5(^), 5(^), where the scale k is comparable to Aasp on the curve a = aoo{S) dehne 
critical covariant gauges (CCG). Remarkably, we only hnd real intersections 
in the asymptotically free regime when the number of quark flavors does not exceed the 
number of colors 

Hf < n (4.8) 


Note that for nf > n, < /3o for any real value of 5. The scale k in (O) is therefore 

irrelevant compared to Aasp for > Uf > n. In the space of parameters 6 ,a,nf/n the 
condition ( [4.6|) is valid on a surface which intersects the surface where ( p.3|) is satished on 
a critical curve. Projections of this curve onto the {6, a)- and (n//n, 5)-planes are shown 
in Fig. 1 . Only on this critical curve in the parameter space, k and Aasp are both asymp¬ 
totically relevant scales and proportional to each other. The perturbatively inaccessible 
region is shaded grey in Fig. 1 and given by na^oS'^ > Po, since k S> Aasp in this case. We 
emphasize that (|4.8|) and (|4.7]) are consistency requirements for a perturbative evaluation. 
In other gauges (or when rif > n) there is no non-trivial hxed point of the moduli 7 . In 
the grey region of Fig. 1 perturbation theory is not consistent and non-perturbative held 
conhgurations are perhaps important even asymptotically. In the white region of Fig. 1, 
the scale k is irrelevant and the only hxed point of our moduli-space is the trivial one. 
The perturbatively interesting region is the phase boundary where k ~ Aasp- 




Fig. 1: Projection of the critical curve onto the (5, a)- and {5,nf/n)- planes. The two points 
on this critical curve for a fixed value of rif/n < 1 define CCG. The perturbatively inaccessible 
region in the parameter space is shaded grey in the (<5, nj/n)-plane. The origin of the axes is at 
the critical point {nj = n, 6 = 1, a = 3). 


The gauge group of QCD is SU{3) and there are rif = 6 quark havors. Since u/ = 6 > 
3 = n in this case, the perturbative analysis indicates that k is an irrelevant asymptotic 
scale for QCD. However, only 2 — 3 light quark havors seem to be important for the low- 
energy dynamics of QCD and heavy quark havors are furthermore unstable when havor- 
changing weak interactions are included^. It is not inconsistent phenomenologically that 
K is irrelevant at energies much above the J/T threshold where more than 3 quark havors 
are dynamical and where ordinary perturbation theory works reasonably well. 

^^The coupling to photons invalidates our perturbative analysis, which depends heavily on asymptotic 
freedom 
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To qualitatively describe the dynamics of QCD below a few Geld, an SU{3) gauge 
theory with only 3 or even less dynamical quark flavors appears to be a reasonable model, 
especially when weak interactions are ignored. Heavy quark flavours appear to be ade¬ 
quately described by heavy quark expansions |T 8 [| around p/ ~ cx). The purely gluonic 
theory with Uf = 0 and/or the large n limit of SU{n) gauge theories are furthermore 
of theoretical interest, since Wilson’s criterium[Q for absolute conhnement is only valid 
in the absence of dynamical quarks. In these limiting cases, k is an asymptotic scale in 
CCG. 


Let us therefore explicitely verify that ( 7 ) is a truly marginal parameter in a (hypo¬ 
thetical) SU{n) gauge theory with nj < n quark flavors in CCG. The critical exponent of 
7 vanishes if the anomalous dimension of 7 ^ is of order The renormalization constant 
Z ^2 relating bare and renormalized moduli, is calculated in Appendix B to 

order g'^. For arbitrary gauge parameters 6 , a one obtains 

z,. = 1 - (3 + a(l - 2^)) + 0(y) (4.9) 


In CCG dehned by the hxed points ( p.3|) and ([4.7|), the coefficient of g"^ in 
vanishes and the critical exponent of 7 ^ thus vanishes as well. We thus have explicitly 
verihed the relation between the two coefficients in ( |4.2|) at hxed points of the parameter 
space. The non-trivial absolute minima of the effective potential dehning the scale k 
therefore do not depend on the value of the coupling g and dehne an asymptotic RG- 
invariant scale k which is a physical mass scale in GGG[^ As far as the perturbative 
evaluation of correlation functions invariant under the global SU (n) symmetry in GGG is 
concerned, 7 in ( | 2 . 11 |) can be diagonalized and its eigenvalues treated as mass-parameters 
proportional to A^gp that do not depend on the coupling. 


As remarked earlier, the absolute minimum of the one-loop effective potential (|4.4|) 
is unique modulo SU{n) transformations. It is attained when none of the eigenvalue 
differences vanish, and the SU{n) group is therefore maximally broken to U^~^. Since 
the difference of two eigenvalues is proportional to the asymptotic scale parameter and 
furthermore RG-invariant, a variant of our previous argument indicates that the absolute 
minimum of the full effective potential also only possesses a invariance in GGG. 

Higher order- and non-perturbative- corrections could only lift degeneracies among the 
eigenvalues and further break the symmetry. The difference of any two eigenvalues (for 
instance the lowest two) at the hxed point of the moduli space however dehnes a physical 
and asymptotic scale in GGG which cannot vanish in the full effective potential either. 


We next verify our previous assumption that ghost number is conserved, i.e. that a 
non-trivial hxed point with {(pa) 7 ^ 0 would be perturbatively inconsistent. We only need 
to show that the critical exponent of a(p does not vanish anywhere on the curve aoo{^) 
of (|3.3|) for Hf in the domain (|3.7|) . A one-loop calculation similar to the previous one 


^•^Note that for the extreme case uf = n, the two fixed points coalesce at Uoo = 3, <5^ = 1, a kind of 
Yennie-gauge[^^ with nice infrared behavior of the quark propagator]^, and vanishing critical exponents 
of the ghost and gluon helds 
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gives for the renormalization constant relating the bare and renormalized helds 

z,* = 1 + M(1 - 2i) - 3) + O(J-) (4.10) 

The critical exponent of ((T0) therefore vanishes only if aoo^(l — 25) = 3. Demanding also 
(|3.3|) , this occurs only when the gauge parameter 5 is a solution of the quadratic equation 

5(1 - 25)(13n - Anf) = 9n(l - 25 + 25^) (4.11) 

( [4.111) which has real solutions only for 


-5 + 18^ 13 

Hf > -:- n > —n 

4 4 


(4.12) 


But Uf > is outside the range (|3.7|) for which a approaches the hxed point aoo(5) 


13 „ 


asymptotically. Furthermore ( |4.10|) is not of order at the hxed point a = 0 for rif > 
either. The assumption that ghost number is broken therefore would be perturbatively 
inconsistent. 


Consider hnally the moduli a^, the constant parts of the gauge held. For a non¬ 
trivial hxed point with an anomalous dimension of order the term of order g'^ of the 
renormalization constant Za = Z^Z"^ would have to vanish. From (|3.2|) we see that this 
only occurs for a = —3 - clearly an impossible value for the gauge parameter. When 
a > 0, the anomalous dimension of is negative, and the only perturbatively stable 
hxed point is = 0. 

To summarize: the perturbative analysis of this section indicates that the global 
S'[/(?7,)-symmetry is spontaneously broken to [/(l)”"^ for Uf < n quark havors in CCG. 
Ghost number on the other hand appears to be conserved. The non-trivial expectation 
value (^ 57 ) is characterized by a RG-invariant scale k at two critical points of the param¬ 
eter space 

(^ = 0,a = aoo(5^),5 = 5^) (4.13) 

given by (|^ ), and (|4(7|) . We next investigate qualitative implications for the low energy 
dynamics of the theory from the spontaneous breakdown of the global SU{n) symmetry 
in GGG. 


5 The broken phase 


Since the equivariantly gauge hxed SU (n) theory is covariant and translationally invariant 
in the inhnite volume limit, spontaneously broken symmetries should give us access to 
the non-perturbative infrared dynamics of the model by Goldstone’s theorem[pl[|. 


A non-trivial expectation value ( 7 ) 7 ^ 0 does not imply a breakdown of the equivariant 
BRST-symmetryp^. (sc(x)) in fact is the conjugate variable to 7 in (|2.11|) in the sense 

^^The question of a global anomaly in the BRST-symmetry of the equivariantly quantized continuum 
model for certain base manifolds was addressed inj|]. 
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that 


(5.1) 


OVettK'l) 

djO- 


= (sc“(a;)/c/)^ = (7“/(ac/^) - 6f 




where we have used the equation of motion for 7 and translational invariance. A nec¬ 
essary condition for an unbroken equivariant BRST-symmetry therefore is that 7 is an 
extremum of 14//. Conversely, if we require an unbroken BRST-symmetry ( |2.13| ), the 
possible expectation values for the moduli 7 are solutions to the gap equation 


/ 7 


(5.2) 


Note that ( |5.2|) in general is a somewhat weaker statement than that ( 7 ) correspond to 
an absolute minimum of the effective potential. In the last section we saw that non-trivial 
RG-invariant solutions to (|5.2|) exist in CCG. 


An superhcially similar equation to ( b. 2 |) has previously been obtained in the thermo¬ 
dynamic limit of lattice gauge theory in Landau gauge as a necessary horizon condition 
for effectively constraining the gauge theory to the fundamental modular region|^. Apart 
from the fact that (|5.2|) also introduces a mass scale via a consistency requirement, this 
gap equation differs from the horizon condition in many respects. Perhaps the most ob¬ 
vious difference is that the Landau-type gauges of@] cannot be investigated in our way, 
since they do not correspond to perturbatively stable hxed points of the gluonic theory. 
Furthermore, a non-trivial solution to (|5.2|) breaks the global SU{n) symmetry, which was 
assumed not to be broken in the lattice gauge theory |p. Finally, it was shown in^] that 
a non-trivial solution of the horizon condition spontaneously breaks the BRST-symmetry 
of the lattice gauge theory in the thermodynamic limit. The gap equation ( |5.2| ) in our 
case on the other hand is a minimal requirement for the equivariant BRST-symmetry of 
the model to remain unbroken. 


Although ( b.2D is only a necessary condition, we assume in the following that the 
equivariant BRST-symmetry is not broken in GGG. Perturbation theory gives no indi¬ 
cation to the contrary and the breaking of BRST found inP| for certain base manifolds 
does not occur for a torus. It should however be noted that the naive BRST-symmetry 
of covariant gauges p 0 | is broken - it would give a vanishing partition function and was 
the reason for restricting to an equivariant cohomology. It is therefore perhaps not so 
surprising that we hnd non-trivial hxed points of the moduli space with this equivariant 
construction in certain gauges. 


(^.2|) relates ( 7 ) to the expectation value of a local composite operator. To hnd 
the Goldstone bosons of the spontaneously broken symmetry, we hrst break the SU{n) 
symmetry explicitly and extend the action (|2.11|) by a BRST-invariant local source term 


So — S[p] = So -\-2a J Tt p{x)sc{x) 


(5.3) 


We now study local Ward-identities of this extended model in the limit of vanishing almost 
constant source p{x) in the inhnite volume limit. 
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To derive the relevant Ward-identity, consider the following local inhnitesimal variation 
of the helds in the functional integral 


SAf,{x) = D^e{x) = d^9{x) + g[Af,{x), 9{x)] 

Sipiix) = -g9{x)'ipi{x) 5il)i{x) = g9{x)'ilji{x) 

Sc(x) = g[c{x), Sc(x) = g[c{x), 9{x)] 

Sj = Sj = 6a = Sa = Scj) = 0 


(6.4) 


The transformation of the fermion- and gauge-helds in ( |5.4|) is just an inhnitesimal gauge 
transformation and thus only the gauge hxing part of the action (|5.3|) changes to hrst 
order in 9{x) 


69°'{x) 


The variation (|5.4|) gives the color current 


= —df^{J^{x) + D'^^p^{x)) -f -|- a6g'^f’^^^c^{x)c^{x))p^{x)- 

—— g5^^ f^^c'^{x)c^{x) + a^c^{x) — 

-l(jhfcde^d(^x)c\x) - f^'^^d^{x)c\x)) 


J<4{x) = 


-sD^^c^{x) — a (x) = 


= a ^{d■ A^{x)) -I- g6f‘^^^&{x)D^d^{x) — g{l — 5)f°'^^A’{x)D"^c'^{x) 


(5.5) 

(5.6) 


without use of the equations of motion. The current ( ^. 6 |) does not explicitly depend on 
the moduli and in fact is just the usual one for gauges with quartic ghost interaction!]^. 
Contrary to the standard case however, this current is not s-exact for 7 7 ^ 0. Note that 
the constant ghosts do not transform under the local variation (|5.4|) and that the change of 
the action is therefore not proportional to the divergence of the color current (|5.6|) . This 
somewhat complicates the Ward identities, but we will see that the additional terms are 
irrelevant at vanishing momentum transfer. Because the gauge symmetry is not anoma¬ 
lous, the measure of the path integral is invariant under the change of variables ( b.4|) . We 
thus have the Ward-identity 


' 50 ' 
^ 69°'{x) 


O 


69°'{x) 


= 0 


(5.7) 


where {■ ■ indicates that the expectation value is for hxed source p{x) and O is some 
functional of the helds. With O = 0°‘{x) 

0^{x) = ga6r^^c\xy{x) (5.8) 

(^) at vanishing momentum transfer in the limit of inhnite space-time volume becomes 


/A (^0%x)d-{J\y) + D’’‘^p%y))) = (^/“"''C>"(x)) + 


+^/A/'’"V(2/) 


V(2/) 

d 


(0“(a;))^ 
d 


(5.9) 
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+ ( O’^ix) gp^^ 7^— -1- 7^=— + a^— + a 
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da^ 
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Since C>“(x) does not explicitly depend on constant ghosts, the last term in O is 
proportional to equations of motion of the moduli and vanishes. In the limit p{x) —0 of 


a slowly varying source, (|5.9|) requires that 

jd*x (d-r{x) 0%)) = jd*x (a-J°(x) gaSS'^if{,y)c%)) = 


(5.10) 


where we used (|5.1| ) to express the expectation value of the operator {0°'{x ))in terms 


of 7 “. Since the RHS of ( |5.1CI|) does not vanish at the non-trivial minimum of the effective 
potential in CCG, (|5.10 ) implies a long range behavior of the correlation function 


j;{x) o\Q) 


f 


abcgx^c 


X,. 


( 6 . 11 ) 


2ti‘^x* 

which is the signature of massless Goldstone modes. Note that the last step assumes 
translational and Euclidean invariance of the thermodynamic limit. 


The relation ( |5.11|) between ( 7 ) and the correlation function of the composite operator 
0°‘{x) of (1^ ) with the color current (jT^) does not appear to be RG-invariant, since 0°‘{x) 
generally mixes with d-A. It is however straightforward to check that the divergences 
on the LHS of (|5.11|) cancel to order g'^. We have observed earlier that = 1 + 0{g‘^) 
precisely at the hxed points ( 4.13 ) dehning GGG. The relation ( 5.11 ) is therefore consistent 
in the weak coupling regime of GGG. 


However, divergences in the correlation function (|5.11|) only cancel to hrst order in the 
loop expansion for the combination of composite operators in (b. 6 |) . The mixing between 
individual terms in implies that the Goldstone pole must also be present in every 
correlation function of 0°'{x) with any composite operator appearing in the color current 
and in particular in the IPR part of 


d^d-A'^{x) O\0) 


2 x^Z(/i) . 

-J 7 


x^ 


(5.12) 


The residue Z (/i) in (|5.12|) depends on the renormalization point p even at arbitrary weak 
coupling - only the residue in ( 5.11 ) of the Goldstone contribution to the correlation of 0°‘ 
with the full current is RG-invariant in GGG in the limit p —>■ cx). Z{p) however cannot 
vanish for arbitrary values of the renormalization scale, because the operator d^d-A°' is 
necessary for a cancelation of the one-loop divergences on the LHS of ( |5.11| ). A coupling 
to the Goldstone modes is therefore inevitably generated in (|5.12|) by a change of the 
renormalization scale. 


In the translationally invariant theory, ( 5.12|) requires a long range behavior of the 
■‘mixed” propagator of the form 


d-A’^ix) O\0) 


■^(h) jabc-c 


X^ 


(5.13) 


We now argue that a long-range behavior of the longitudinal gluon propagator 
'd-A%x) a-A'’(0)) ^ ^Mjacdjbce-d-e 


X^ 
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and of the C>“-propagator 


(0%X) C>'’(0)) 

also is a natural consequence of the massless Goldstone excitations. In the case of covari- 
antly quantized gauge theories this coupling to Goldstone excitations is unfortunately not 
as inevitable as in a theory with a positive dehnite Hilbert space. 


The residues in (|5.14|) and (|5.15|) are related by the s-symmetry (|2.13|), which implies 


that 


(sc°‘{x)sc^{y)'^ = — (^c°'{x) ss(^{y)^ = a — y) 


(5.16) 


In gauges <5 = 0 this is the Slavnov-Taylor identity which determines the longitudinal 
part of the gluon propagator completely. In general gauges with <5 7 ^ 0, (|5.16|) only 

relates connected Green’s functions and in particular requires the absence of Goldstone 
excitations in the combination 


0 = a^ 


(O'^ix) - (C>“)) (C>'’(0) - {O'’))) + {d-A'’{x) d-A'’{0) 


+ {{0^{x) - (C>“)) a-H'’(0)) + (d-A%x) (C>'’(0) - 


x ^>0 


(5.17) 


for x"^ > 0 . 


The color structure of (|5.13| ) leads to a cancelation of the long range parts of the 
last two terms in (|5.17|) . The long range correlations of the hrst two terms also have to 
cancel for ( b.lTD to hold. Note that for S ^ 0, (^.17|) does not demand vanishing residue 
Zl, because the metric of the Hilbert space of the covariant gauge theory is not positive 
dehnite. 


States of negative norm in covariant gauges however also prevent us from proving in the 
usual way that in ( 5.14 ) does not vanish. A non-trivial cancelation of the Goldstone 
poles between the hrst two terms of (|5.17| ) however only requires that the vertex 
between the scalar cc-mode and longitudinal gluons with an antisymmetric color structure 
does not vanish for —> 0. We perturbatively verify that this is the case in the broken 

phase of the S'[/(2)-model in Appendix G. The one-loop calculation of this vertex gives 


p[a6] 

oaA 


(p^ —>• 0) = £ 


abc-c 


2 x2 


2>2tik^ 


(6.18) 


Note that this dimensionless antisymmetric coupling arises only when the global SU{2) 
symmetry is spontaneously broken. A vanishing of the full O vertex at zero momentum 
and weak coupling would require a cancelation of this leading contribution from non- 
perturbative corrections that are of order 1/p^ cx ln(A^/p^). Such a dependence on the 
renormalization scale /i is generally induced only by divergences. Since the vertex (|5.18|) 
is not present in the action ( |2.11|) and the model is renormalizable, the (9 cl ■ A-vertex 
antisymmetric in color is hnite to all orders. We thus do not expect a cancelation of the 
lowest order vertex (|5.18D at weak coupling. 
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The color structure of ( p.l8|) guarantees the cancelation of Goldstone poles in the 
Ward identity (|5.17|) . Since the adjoint representation is the only one common to the 
antisymmetric tensor product of two adjoint representations and the symmetric tensor 
product of any number of adjoint representations, the corresponding vertex for SU{n) 
can only be proportional to , where + p 2 d '^'^^+ ... 

is an odd and totally symmetric polynomial in 7 in the adjoint representation of the 
group. That the color structure of the Goldstone contribution to {d-A[x) 9-y4(0)) is of 
the form (|5.14| ) then follows from ( b.llj| ) and Bose-symmetry. The real, symmetric and 
positive semidehnite — 1 ) x — 1 ) matrix exactly (n — 1 ) 

vanishing eigenvalues when 7 breaks the S't/(n)-group to 17(1)"'“^. Thus the n{n — 1) 
positive eigenvalues of and the number of expected Goldstone modes correspond. 
Of some importance for the next section is that the vertex Toqa is of order for large 
n, since the expectation value ( 7 “) is of order and insertions of ( 7 ) therefore do not 
change the leading order in n of a diagram. The one-loop diagram shown in Fig. 5, which 
we evaluated for SU{2) in Appendix G, is planar and thus of order in SU{n). 


Finally note that the critical exponent of the residue Zi{fi) in (|5.14|) is given by the 
anomalous dimension of the gluon held of order since the anomalous dimension of 7 
is of order g'^ in GGG. 


Although indicative, we are well aware that these arguments for a Goldstone contribu¬ 
tion to the longitudinal gluon propagator in GGG are not conclusive. In the next section 
we therefore present further evidence for a long range behavior of the longitudinal gluon 


propagator of the form (|5.14|) : assuming that the dominant long-range contribution to a 
Wilson loop in GGG comes from this Goldstone contribution to the longitudinal gluon 
interaction, the expectation value of large Wilson loops is found to decay exponentially 
with the enclosed area in the planar gluonic theory. 


6 The Wilson loop in CCG 


The order parameter of conhnement proposed by Wilson [^] is the expectation value of the 
non-local functional W (C) 


W{C) := TiV exp 



( 6 . 1 ) 


where V denotes path ordering and C is a generic smooth path. In a conhning gluonic 
theory [nj = 0) the expectation value of ( | 6 . 1 |) should decay exponentially with the area 
of the minimal surface enclosed by C = dD 


In {W{dD)} - aArea(D) (6.2) 

where the constant of proportionality a is the string tension. This should in particular 
be the case for discs D with large Euclidean radius R. 
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{W{dD)) is a gauge- as well as RG-invariant observable for any SU{n) gauge theory. 
It suffices to show that the string tension a does not vanish for large n, i.e. for the subset 
of planar diagramsin the expansion of the Wilson loop. A cancelation between planar 
and non-planar contributions to the expectation value precisely for n = 3, the gauge group 
of interest, would be unfortunate. We will therefore only study planar contributions to 
the expectation value (|6.2|) of circular Wilson loops such as the one shown in Fig. 2. 



Fig. 2: A planar gluonic contribution to a circular Wilson loop. 


Although the leading order in n of (W{dD)) is gauge invariant, the contributions 
from planar connected gluonic Green’s functions to a diagram such as Fig. 2 are gauge 
dependent individually. It is therefore possible, that some of these contributions to the 
expectation value are negligible compared to others in certain gauges when the area of the 
loop is large. In the preceding sections we argued that the longitudinal gluon propagator 
in GGG is of particularly long range due to Goldstone modes0. The planar connected 
Green’s functions contributing to large Wilson loops naturally fall into two categories in 
GGG: i) those which are dominated by the Goldstone poles at long range and ii) those 
which aren’t. Gontributions of type i) at large distances can be regarded as longitudinal 
gluon exchanges with an effective coupling G(/i) to the loop. G(/i) coincides with the 
gauge coupling g{fi) in lowest order of the weak coupling expansion only and is also of 
order 1/s/n. (The order g^{yi) corrections to G{g) are obtained in Appendix D. One 
hnds that G(/i), the effective coupling of a longitudinal gluon to the Wilson loop, and the 
perturbative coupling g{g) differ in their dependence on the renormalization scale /r.) In 
critical gauges contributions of type i) give rise to a long range interaction 




oo) = G\g) lAlix) A^(0) 


Goldstone 


(6.3) 


from the exchange of n{n — 1) Goldstone bosons. Because the SU{n) group is broken to 
U (1)"^“^, the number of Goldstone bosons is of order n^. To show that the interaction (|6.3|) 
is also of leading order in n, however requires that their average residue does not vanish 
in the limit of large n. As remarked earlier, it is evident from the action ( 2 . 11 ) that 
insertions of 7 are of order and thus do not change the leading order in n of a diagram. 
Thus the effective potential (|4.4|) is of order n, but its minimum, ( 7 ) as well as the scale 
K are of order n^. Since all eigenvalues of ( 7 ) differ and each is of order n°, the sum of 
their squares therefore grows like n, that is 


^^With n/ 


^a-a ^ _ 2 Tr 72 = = 0 (n) 

i=l 


0 the two CCG are at gauge parameters <5^ = q,± _ 


(6.4) 
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because at most one of the eigenvalues vanishes. To leading order in n we can therefore 
replace the color dependence of (|5.14|) by the average 




-- 

n 


+ 0{l/n) 


(6.5) 


Due to ( |6.4|) , the Goldstone contribution to the longitudinal propagator is of 0(1), if the 
residue Zioi the Goldstone modes is of order n°. The previous argument that the coupling 
Zl does not vanish is also valid in leading order of n (since the Ward identity ( |5.10| ) is 
of leading order and the leading contributions to the vertex T^^^ shown in Fig. 5 are 
planar). As far as color counting is concerned, the exchange of Goldstone modes in GGG 
therefore contributes to the Wilson loop in leading order of the 1/n expansion. 


Integrating ( |6.5|) we obtain the long-range behavior of the longitudinal gluon propa¬ 
gator to 0(1) 

Ki^)) - ZL{fi)6^^^^d^d^x‘^ln{b{n)x^) ( 6 . 6 ) 

where we have assumed the least singular behavior for —>■ oo compatible with (|6.5|) . 
The integration constant 6(/i) cannot be determined at this stage. We will numerically 
hnd that 6(/i) is related to ZlG^ if the physical spectrum is unitary. At present we only 
remark that changing the parameter b alters the propagator by a constant and therefore 
changes its long range behavior. The Goldstone pole in the correlator ( |6.5|) , which was 
inferred from the Ward identities, only determines the leading long-range behavior of the 
propagator which does not depend on the associated scale b. Note however that an explicit 
solution of the model, would give the long range part of the gluon propagator up to a 
physical scale and determine the associated dimensionless parameter Zi/b uniquely in 
terms of the dimensionless coupling constants. 


Inserting (1^) in (i]D we arrive at a parameterization of the long range 2-point 
interaction in leading order of n 


,ab t 


^ab 


vZJx^ ~ oo) = — K^dnd„x^ \n(sK^x^/e^) 
^ n 


(6.7) 


in terms of K and s which are related to G^Zl respectively b/{G‘^ZL). In (|6.7|) s is a 
dimensionless parameter and K is an inverse length. In Appendix D we verify that these 
parameters approach renormalization group invariants for /i —cx). To establish this, we 
consider the perturbative planar 1-loop corrections of Fig. 6 to the effective long-range ex¬ 
change ( |6.7|) within a Wilson loop. The explicitly calculation of Appendix D shows that 
the UV-divergences associated with an evaluation of diagrams 6a-|-b precisely cancel the 
renormalization factor for the propagator and coupling to order g^. No UV-subtractions 
of order g‘^ are therefore necessary to render the contribution to the expectation value of 
the Wilson loop from effective planar longitudinal 2-point exchanges hnite. The interac¬ 
tion (|6.7|) therefore does not depend on the renormalization scale p as p —> cx). 


In Appendix E we prove that planar contributions to the Wilson loop generated by an 
effective 2-point longitudinal interaction v are resummed by a nonlinear integral equation. 
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In the case of a circular Wilson loop of radius R this integral equation is of the relatively 
simple form 

Wn{9) = 1 + r dO' d9’'vR{9'')WR{f)WR{9' - 9") (6.8) 

JO JO 

In (|6.8|) the color traces have been performed and the expectation value of the circular 
Wilson loop W{BDr) of large radius R with the above approximations is 


{W{dDR)) nWR{2'K) . 

The effective interaction vr{9) in (^) is related to by 


(6.9) 




89X 89y n 
1 8^ 


2 89 ‘^ 


+ 0(l/n^ 


{KRf sm\9/2) \n{s{KR/eY sm\9/2)) 


= {KRy 1 + cos(6') ln(s(i^i?)^ sin^(6'/2)) 


( 6 . 10 ) 


since the distance {x—yY = AR^ siv?{{9x—9y)/2) of two points on the loop depends only on 
the radius R and the angular difference 9 = 9x — 9y between them. The interaction ( 6.1CI| ) 
depends parametrically on the dimensionless area {KRY of the loop and the dimensionless 
parameter s. 


We numerically solved the nonlinear integro-differential equation 


^Wr{ 9) = f d9'vR{9')WR{9')WR{9 - 9 '), 
o9 JO 


IVr(O) = 1 


( 6 . 11 ) 


which is equivalent to ( |6.8|) . For small values of 0 <C {KR^ (|6.11|) was solved iter¬ 
atively and we then used these starting values to extend the solution to 0 = 27r by a 
predictor-corrector algorithm based on the Adams-Bashforth-Moulton scheme. The nu¬ 
merical integrations where performed with a modihed Rhomberg algorithm||23||. We found 


the predictor-corrector method appropriate in this case, since the integral in ( |6.11 ) de¬ 
pends on all previous values for the function Wr. Wr{9) itself is a reasonably smooth 
function which oscillates a few times in the interval [0,27r]. Special care was taken to 
correctly take into account the logarithmic (but integrable) singularities oi vr aX 9 = Q 
and 9 = 2tt. The algorithm is exceptionally stable and we quote results which we believe^ 
to be accurate to ~ 10“^^. The results of these numerical simulations are summarized in 
Fig. 3. For a particular value s = Syrit. = 1.195, we numerically hnd that the expectation 
value of the Wilson loop decreases exponentially with the dimensionless area {KR)'^ over 
12 orders of magnitude. The results for the Wilson loop with neighboring values s = 1.17 
and s = 1.25 are also shown in Fig. 3 for comparison. We observe that for s = 1.17 < Scrit. 
the logarithm of the Wilson loop decreases faster than linear with the area. Such a be¬ 
havior of the expectation value is excluded by causality 
also observe a violation of causality 


For s = 1.25 > Scrit. we 
for large values of the area, the expectation value 


used up to 10^ integration points in the interval [0,27r] and the predictor-corrector algorithm is 
of fourth order. 
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increases again. This behavior is more marked at even larger values of s, as can be seen 
for s = 2.0 in Fig. 3. At s = Scrit. the oscillations which eventually cause the expectation 
value of the Wilson loop to become negative are of order 10“^^ and thus well below the 
numerical accuracy of the algorithm. 


The results summarized in Fig. 3 conhrm our expectation that the long range inter¬ 
action (|6.7|) is only consistent with unitarity at a particular value of the dimensionless 
parameter s. Our best numerical determination of this critical value for the parameter 
is s = Scrit. ~ 1.195 ± .01. A full solution of the model would predict a value for the 
dimensionless parameter s in the asymptotic interaction and furthermore relate the scale 
K to A ASP- With the assumption that the approximations made in the evaluation of the 
expectation value of large Wilson loops are reasonable in CCG, the numerical evidence 
suggests that the string tension does not vanish and gives an estimate for Scrit. if we re¬ 
quire unitarity. Note that Scrit. is a dimensionless parameter which would be a prediction 
of the gauge theory, whereas K is & scale, which essentially has to be determined experi¬ 
mentally. If the residue in (|5.14|) and the expectation value 7 can be related to A asp 
in a two-loop calculation, a quantitative estimate of Aasp/K eventually could be possible 
within the framework of perturbation theory. 




Area 


Area 




Area 


Area 


Fig. 3: The expectation value of large circular Wilson loops as a function of the enclosed dimen¬ 
sionless area {KKf‘ in the long range approximation of CCG. 

Top-left: The logarithm of the expectation value for s = 1.195 in the range {KR)‘^ < 4; the 
expectation value is positive in this range. The line at log((IF)) = —12 is an estimate of the 


27 









accuracy of the numerical solution. 

Top-right: same as top-left but for the greater range (KR)^ < 8; the plot shows the logarithm 
of the absolute value of the expectation value, which changes sign at the sharp dentures for 
(KR)^ > 4. We attribute these oscillations of order 10“^^ in the expectation value to numerical 
noise of the algorithm. 

Bottom-left: The logarithm of the absolute value of the expectation value of the Wilson loop in 
the range {KR)‘^ < 4 for parameter s = 1.195 (bold curve) as well as s = 1.17 and s = 1.25; for 
s = 1.17 the logarithm of the expectation value decreases faster than linearly with the area and 
violates unitarity. 

Bottom-right: Same as bottom-left for a wider range of the parameter s and for {KR)‘^ < 8. The 
result for s = 2.0 shows a violation of unitarity at the level of 10“^ which could be interpreted 
as due to a low-mass state of negative norm. 


Our numerical studies also indicate that the Goldstone pole in ( 5.14 ) is essential for an 
exponential decrease of the Wilson loop in this approximation. The logarithmic depen¬ 
dence on the distance of the effective interaction (|6.7|) and ( |6.1(]| ) is a direct consequence 
of, and implies, a Goldstone pole in (|5.14|) . We numerically investigated the effect of 
replacing the logarithm in ( p.l0| ) by a constant and found that this immediately leads to 
a non-unitary dependence of the expectation value of the Wilson loop. Our numerical 
results therefore depend sensitively on the asymptotic effective interaction ( p.lOj ) and are 
not a generic consequence of the non-linear integral equation ( |6.8|) . Partly due to this sen¬ 
sitivity however, we presently do not analytically understand the solution to this nonlinear 
integro-differential equation. An analytic, or approximate, solution to (|6.11|) for large val¬ 
ues of {KRY would provide a better understanding of this conhnement mechanism in 
GGG. 


7 Discussion 


Since quarks have not been observed as asymptotic states, conhnement is one of the 
central issues in SU{n) gauge theories without fundamental scalars. On the other hand, 
these models are asymptotically free and perturbation theory in conjunction with the 
RG-equation is a powerful tool in their investigation. It was only recently observed that 
a translationally invariant BRST-quantization of these theories on compact space-time 
manifolds introduces a non-trivial moduli space |^, Q because the BRST-symmetry has to 
be realized equivariantly with respect to global gauge transformations. In this work we 
found non-trivial hxed points on this moduli space in GGG, which respect the equivariant 
BRST-symmetry but break the global S't/(n)-symmetry spontaneously to [/(1)"“^. The 
associated Goldstone bosons are cc-bound states in the adjoint representation, which we 
argued lead to an IR-singular long distance behavior of the longitudinal gluon correlation 
of the form ( |5.14|) . Our numerical results for the Wilson loop indicate that this long range 
behavior of the longitudinal gluon propagator conhnes. 


This scenario is however only realized in particular covariant gauges and only for np < 














n quark flavors. In gauges other than CCG, the behavior of the Wilson loop is presumably 
similar, but the dynamics leading to conhnement cannot be described in this way. More 
specihcally, the relevant unphysical degrees of freedom that give a conhning behavior 
of large Wilson loops in other covariant gauges are held conhgurations which cannot 
be described by the limited moduli-space at our disposal. The situation is somewhat 
analogous to atoms in QED: the spectrum of an atom is gauge invariant, but its qualitative 
description in gauges other than Coulomb gauge is a non-trivial matter. In Coulomb 
gauge, the unphysical degrees of freedom conspire to give a potential which describes the 
spectrum qualitatively and radiative corrections are computable and small. We believe 
that gauge freedom can also be exploited advantageously in SU (n), but with the important 
difference that these non-abelian models are asymptotically free. Perturbation theory is 
valid only at high energies and we rely on spontaneous symmetry breaking to relate 
perturbative information to low energy dynamics. Within our limited moduli space, we 
can verify such a symmetry breakdown only in CCC and only for Uf < n. As far as we 
know, the expectation value of some composite scalar not described by our moduli space, 
would show that SU{n) is broken to also in covariant gauges other than CCC, 

and the corresponding Coldstone excitations could similarly lead to conhnement. The 
limitation of our perturbative approach is that we cannot access the relevant composite 
scalars in gauges other than CCC. In non-covariant gauges the conhnement mechanism is 
probably quite diherent and we have no suggestion how to perturbatively gain information 
about the low energy dynamics of the theory in this case. 

The outstanding (and new) feature of the proposed conhnement mechanism in CCC 
is the Coldstone pole in {d-A{x) cI-A(O)), which implies that the longitudinal gluon inter¬ 
action is proportional to ~ 1/k* for small momentum transfers. In conventional covariant 
gauges, the BRST-symmetry constrains the longitudinal propagator to be proportional 
to 1/fc^; a l/k'^ behavior of the longitudinal propagator however does not violate Ward- 
identities in the extended covariant gauges with quartic ghost interaction (5 7 ^ 0). We ar¬ 
gued that perturbation theory is consistent with a dynamical spontaneous breaking of the 
global S'17(n)-symmetry for certain critical gauges due to ghost — antighost-coudeusution 
in the adjoint. Remarkably, this spontaneous symmetry breaking of the gauge hxed the¬ 
ory does not give rise to a mass for the longitudinal degrees of freedom of the gauge held 
via the Higgs mechanism and the vertex (^.181) we obtain for the coupling between these 
modes does not vanish for —>■ 0. This is probably because the coupling of the longitu¬ 
dinal gauge held to the massless (composite) cc-Coldstone modes is not constrained by 
gauge symmetry. 

We presented numerical evidence that the observed behavior of the longitudinal propa¬ 
gator in non-abelian theories would conhne in the sense of Wilson|Q. An appealing feature 
of this scenario is that the IR-singular behavior of the longitudinal exchange does not lead 
to a strong “Van der Waals” force between two distinct Wilson loops at large distances, 
because the interaction is a total derivative in both loop parameters. For sufficiently small 
loops, for which all exchanges within the loops can be neglected, the force between two 
Wilson loops due to longitudinal exchanges is therefore proportional to the physical size[^ 

^^Since the contribution from the exchange of one and two longitudinal gluons vanishes, the force is at 
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of either loop. The Van der Waals force from longitudinal exchange between local physical 
observables, for which this approximate evaluation of the long range interaction is valid 
in CCG, is thus seen to vanish. 


Let us finally point out that the symmetry breaking we have discussed gives quantita¬ 
tive corrections to conventional perturbative results. These can in principle be calculated 
and experimentally verified. Of particular interest are perhaps power corrections to phys¬ 
ical correlation functions at large Euclidean momentum transfers. These power correc¬ 
tions are described by vacuum expectation values of composite operators in the operator 
product expansion of Wilsonp^. At present these expectation values are determined 


experimentally or are obtained from lattice simulations. In CCG such power corrections 
arise naturally in the loop expansion due to the expectation values of the moduli. 
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A The partition function (1) 


We here adapt the proof|^ that the partition function of the TQFT of the gauge group 
does not vanish to our case. Consider the partition function 

Z[A] = (1)^ = J[dU][dc][dc][db]d(j)dadad^d'y (A.l) 

of the equivariant TQFT on the gauge group with action 

= sWgf[A^] (A.2) 


where Wgf[A^] is the gauge fixing functional (p.9|) with A replaced by 


A,, 


= U^AM + -U^dM 


(A.3) 


A^ in (|AA|) is a background connection that identifies a particular orbit. U{x) G SU{n) is 


a local gauge transformation and the action of the BRST-operator s on the fields is defined 
as in (p.l|) with the replacement (|A-3|) . The action of s on the background connection A 
and the gauge transformation U individually is given by 


sA^(x) = 0 

least proportional to the sixth power of the physical size of either Wilson loop 
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sU{x) = U{x){gc{x) + uj) 


(A.4) 


s thus effects an infinitesimal variation of the gauge group element U. It is straightforward 
to show that s is nilpotent (s^ = 0). Because the action (|A.2|) is BRST-exact, (|A.1|) is 
the partition function of an equivariant TQFT of Witten type|^ on the gauge group. 
Standard arguments of TQFT show that o does not depend on the parameters a and 
S in (p^.9|) and that ([A.l|) is proportional to the generalized Euler character of the space 


£[A] := {U : U{x) G SU{n),d-A^ = Q}/SU{n) 


(A.5) 


of solutions to the Landau gauge condition d-A = Q modulo global (right-hand) gauge 
transformat ions. 


On the other hand the partition function (1) of the equivariantly gauge fixed SU{n) 
theory is proportional to 

(A.6) 


(1) = j [Vi,Vi,VA]ZlA], 


Ac 


where Sc is the classical gauge invariant action (|2.6|) . Since the gauge symmetry is not 
anomalous, the change of variables 


A 


u 


A 


Uifj —> xjj , ^ '0 


(A.7) 


decouples the [/-integration which is assumed to be normalized. The remaining path 
integral is just the one we are considering and the BRST-symmetry of the TQFT becomes 
the BRST-symmetry (|2.1|) of the model. 

We show that Z[A\ is a non-vanishing constant that does not depend on the back¬ 
ground connection A. Since we are dehning the model on a hnite torus, an explicit mode 
expansion is possible. The Landau gauge condition is associated with fixed points of the 
Morse potential 0 Va[[/] G TZ+ < oo 

Va[U] = [ TtA^-A^ (A.8) 

J M 


which is invariant with respect to global right-hand SU{n) transformations. The gen¬ 
eralized Euler characteristic is the (signed) sum of the Euler characters y of the hnite 
dimensional hxed point manifolds of (|A.8| ) modulo SU{n) 

X(S[A])= y ±x(MilSU(n)) (A.9) 


where the sign depends on the number of negative eigenvalues of the Hessian at a point 
of the hxed point manifold|^. It is clear that this generalized Euler characteristic is 
constant under continuous deformations of the Morse potential Va- Since the space of 
SU{n) connections on a (hnite) torus is connected, Z[A] is thus a constant independent 
of A and it suffices to show that Z[A = A\ does not vanish. To simplify matters, we only 
show that the generalized Euler characteristic of £[A = 0] does not vanish for an SU{2) 
gauge group on a torus. 
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The Morse potential Va=o for A = 0 is also symmetric with respect to left multiplica¬ 
tion by a (constant) group element: U{x) giU^x). Consequently, if U{x) is a critical 
point of the potential (|A.8|) a.t A = 0, 


d-(U^dU) = 0 


(A. 10) 


then so is 

U{x,gL,gR) = gLU{x)gR, '^gx, 9 r^ SU{n) (A.ll) 

Some of these solutions however belong to the same equivalence class modulo right mul¬ 
tiplication by global SUR{n) transformations. We obtain the moduli space of these (right 
equivalence classes of) solutions to ( |A.1(J| ) by noting that left multiplication of U{x) by 
G SU{n) is the same as right multiplication by 

gR{x) = U\x)gLU{x) (A.12) 

giU^x) therefore belongs to the equivalence class modulo right multiplication of U{x) by 
constant group elements only if dgR^x) = 0, i.e. 

[U{x)dU\x), gx] = 0 (A.13) 

Thus left multiplication of U{x) by any gi G SU{n) belonging to the subgroup which 
commutes with U{x)dW{x) is redundant. For an 517(2) gauge group there are only three 
possible subgroups to consider: 


1) gx e SUl{2) satisfy (|AT3D ^ x{SU l{2)/SUl{2)) = 1 

2) g^ G 17(1) C 517z.(2) satisfy (|ATl ^ x{SUl{2)/U{1) - S 2 ) = 2 

3) gx G {1, -1} C SUx{2) satisfy (TO) ^ x{SUx{2)/{l, -1} ~ 50(3)) = 0 

(A.14) 

Note that case 1) implies that U{x)dW{x) = 0 and therefore corresponds to the equiva¬ 
lence class of the identity (modulo right multiplication). There is only one such class. For 
an 517(2) gauge group, case 2) only occurs if U{x)dW{x) = d9{x) is an abelian connec¬ 
tion. Since U{x) should furthermore satisfy (|A.10|) , we can conclude that this connection 
maps the 1-cycles of the torus to 17(1). A typical transformation of this kind on the torus 
has the form 

U{x) = , with [a^, af\ = 0 (A.15) 

with a constant abelian connection a^, whose components on a symmetric torus of exten¬ 
sion L are 

= 27iid-dk^/L, & 7 a (A. 16) 

where the integers k^ are the winding numbers of the mapping of the torus to the 17(1) 
subgroup and d is a unit vector which is the same for all components. Note that the 
corresponding constant abelian connections 

A^{x) = -U^df^U = -a^ (A.17) 

9 9 
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are the Gribov copies of A = 0 on a torus we discussed in section 2.1. They come in 
pairs (a, —a) that are related by a discrete Z 2 subgroup of the S'17ij(2)-transformations. 
A direct consequence is that the expectation value (a^) vanishes for the equivariantly 
gauge hxed SU ( 2 ) theory on a torus. 


From ( |A.14|) we obtain that the moduli space of solutions to ( |A.10|) of an SU (2) theory 
has the topological structure 


Sa=o = SUr{2) X 1 + A 2 X + SO{3) X 


(A.18) 


In other words £a=o/SU (2) can essentially be described as a single point and a collection of 
two- and three-dimensional spheres []^, |^. Although we do not know the Euler numbers 
associated with the topological spaces T and JF, the structure (|A-18|) , together with 
x(S' 2 ) = 2 and x{SO{3)) = 0 suffices to obtain that on a torus 


x{£[A = 0]/SU{2)) = oddj^0 


(A. 19) 


In the vicinity of A = 0 the degeneracy with respect to left group multiplication is in 
general lifted, but the signed sum of Morse indices over the (then isolated) hxed points 
still gives the Euler number ( A.19 ). Since the space of SU{2) connections on a torus is 
connected, we see that the partition function ( |A.1| ) of the TQFT is a constant indepen¬ 
dent of A which does not vanish. To see that the partition function (1) of the gauge 
hxed theory ( |A. 6 |) can then also be normalized, we observe that the determinant of the 
Euclidean Dirac-operator is positive for non-vanishing quark masses in this vector-like 
theory. 


B The one loop effective potential of SU{n) 


We here hnd the dependence of the 1-loop ehective potential on the moduli 7 by extending 
the derivation of the ehective potential]^ for an SU{2) group to SU (n). For reasons given 
in section 4, all moduli except 7 have vanishing expectation values. We take space-time 
to be a symmetrical torus LxLx...xLofD dimensions^. The classical contribution 
to Veffix) is the quadratic term —Tr^'^/ag^ in the action (|2.11|) . The ^f-independent 
term of next order in the loop expansion is obtained by evaluating the inhnite sum of 
one ghost loop diagrams shown in Fig. 4. Crucial for the following is the absence of a 
compensating gluonic one-loop contribution to the ehective potential, because the direct 
coupling Jq- -A'^d-A of 7 to the longitudinal gluon held in ( p.ll|) vanishes on a torus with 
periodic boundary conditions for the gauge held. 

Apart from the finite volume we essentially follow the procedure of ref, using dimensional regu¬ 
larization and the MS-scheme 
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Legend: 

._J.. y- insertion 

ghost-propagation 

Fig. 4: Nonvanishing one-loop contribution to the effective moduli action; a ghost loop with 2p 
insertions of 7 -moduli. 

Note that the antisymmetry of the structure constants implies that only loops 
with an even number of 7 -insertions give a non-vanishing contribution to the effective 
potential. The loop with 2p insertions of the ghost 7 is proportional to 

Cp = TTadj.f" (B.l) 

where we have used the notation X to denote the (anti-hermitian) matrix of the adjoint 
representation of X, Xab = 


2p 


7 r 




V X 


.A 


We want to express the trace (p.l|) in terms of the eigenvalues of 7 in the fundamental 
representation. The trace of interest is given by 


Tr adj.7 = -2Tr 




(B.2) 


k times 


where are the anti-hermitian generators of the fundamental representation of SU{n) 
normalized to Trf“f^ = —Using the Baker-Hausdorff formula the right hand side of 
( [B.2|) is written 


d’^ 

da^ 


Tr 


(B.3) 


Q = 0 


We choose the basis in which 7 is diagonal: ■jkj = i^kdkj (no summation over k). The 
completeness relation for the f“’s, = —2 {dimdij — gives: 

Tr = 1 - iTr e-^“Tr e^“ 

2 2 

The evaluation of the derivatives of this expression for p > 0 leads to the desired expression 
of Cp in terms of the eigenvalues Aj 


2p 


d2p _ _ 2p 

Tr ^ C' 2 ^p(-l)''Tr 7 ^Tr 

k=0 


i:cjy-ir'=i:A‘yAf-* = 2(-i)» y (Ai-Ai 


,2p 


(B.4) 


/c=0 


l< 2 <j'<n 


The contribution to the effective potential from a single loop with 2p insertions of the 
moduli is 


Mp) = 




{2ny) 


E 

{ni...no}^{0...0} 




(B.5) 
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The sum in ( |B.5| ) extends over all the sets of integers {ni.. .rin} describing the complete 
set of modes with momenta = 27rn^/L of the dynamical ghosts. The contribution from 
the constant modes with = 0, ^ = 1,... ,D, to ([B.5D is however eliminated by the 
integration over the ghosts a and 7 in ( |2.11|) . F£,{p) is therefore finite for p > D/A. The 
overall sign of (p.5|) is due to ghost statistics. 


One can analytically continue to non-integer dimensions by casting Fo{p) in integral 
form, 


L^Fnip) = T:T, 




2p 


dxx 


^27r)2y r(2p+l)yo 

where the function f{x) for a; > 0 is the convergent sum, 


2p-l 


-1 


(B. 6 ) 


= Y. 


(B.7) 


Although there is no analytic expression for f{x) at arbitrary values of its argument, the 

^2 


reflection formula 29 


/(^) = V- /(-) 

\ X X 


(B.8) 


implies the asymptotic behavior. 


fix ^ 00 ) = 1 + 2 e-^ + 0 (e-^^) 

2 


/(x^O) = J-{l + 0{e - )) 
V X 


(B.9) 


Since the eigenvalues of the commuting ghost 7 can be treated as real numbers 


, ^ (A, - \i)V 


(47r)2/i4 


(B.IO) 


is positive and the 1 -loop contributions to the effective potential can be summed over all 
values of p > 1. In ([B.IOI ) we have introduced the finite scale p because we are interested 
in the pL —> 00 limit and the eigenvalues Aj in this limit are measured in terms of the 
finite renormalization scale p. Summing over all one loop contributions and adding the 
tree-level potential = Z]i<j( 47 r)^/iSA/na( 52 p 2 gives 




loop 




poo rjrp r . \ D 

/ —sirfi{vijx/2)) Ui'Kx/iLpf)) - 

J 0 X L 


l<2<ji<ri. 


na5‘^g‘^ 


(B.ll) 


The asymptotic behavior of fix) for x —>■ 0 determines the leading behavior of the integral 
in ( B.ll ) when Lp —a 00 . This is the large-volume limit we are interested in. For Lp —a 00 
one thus obtains in H = 4 — 2^ dimensions. 


Vl-loop — p 


(h)V! _ 


cos(-7re)r(£ — 2) 


(B.12) 
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The 1/e-term in the expansion of (|B. 12 |) is canceled by a counter term {1 — Z)L^ 
in the classical action. To order h, we require in the MS-scheme, 


Tr^ 


2, -2e 


Z= 1 + 


nS^ag^jj, 


(B.13) 


Introducing the scale k 


In 


K 


(47r) 


7 b + 1 + 0{\n.g\g^), 


(B.14) 


Aiig? nad^g"^ 

the expression ([B.12|) for the 1-loop effective potential in H = 4 dimensions can be written 






327r2 




(B.15) 


which is ([4.4|) in section 4. 


The counter term of the classical potential together with Za and Zg given in (|3.2|) , 
also determines the renormalization constant Z ^2 of the global ghost 7 to order g‘^. Us¬ 
ing (|B.13|) we obtain in the MS-scheme 


Z;y2 = ZZ^Zg = 1 — 


ng'^jj, 

327r2e 


(3 -|- a(l — 25)) -|- 0(fi'"^) 


(B.16) 


Note that a derivation of the critical exponent of 7 via the effective potential only involves 
UV-divergences and that IR-divergences are absent from this calculation. 


C One-loop contribution to the vertex ( |5.18| ) 


We are interested in the color antisymmetric part of the mixing between the composite 
held 0 °' given in ( |5.3| ) and d-A for vanishing momentum transfer. The one-loop contribu¬ 
tion to the vertex (^.18|) is obtained by evaluating the inhnite sum of ghost loop diagrams 
of Fig. 5 


2 p+l 




0“ y 








Fig. 5: Nonvanishing one-loop contribution to the vertex (5.18). 


Legend: 

y- insertion 
ghost-propagation 


Contributions antisymmetric in the color indices arise only from diagrams with an 
odd number of insertions of the matrix 7 ^^ = For simplicity, we evaluate this 
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contribution to the vertex for an SU{2) gauge group. In a parametric form suitable for 
performing the loop integral, the effective ghost propagator in momentum space is 

1 roo 

Jo 


Gatipy ^ 


p^l + 6'y 

r dXe-^^'Y. 
Jo U 


J'- 


r dXe-^P^ 

-^-h y'f cos(Ak^) -|- 6—— sin(AK^) 

Jo 



where 


-afe _ ^ 


ab x2 7V ^ ^2 laclch 






(C.l) 


(C.2) 


projects on color degrees of freedom transverse to 7 and is the scale of SSB 

of the effective potential {^XA) for an SU{2) theory. The fact that the summation in (G.l) 
gives a relatively simple expression for an SU{2) group, greatly facilitates the evaluation 
of the one-loop contribution to the vertex. Note that the poles of the Euclidean ghost 
propagator (|C. 1 | ) are all at purely imaginary and that the massless ghost is precisely 
that of the unbroken 17(1). The sum of Feynman diagrams with p insertions of 7 in 
space-time dimensions corresponds to the loop-integral 

r?(9) = 

(C.3) 

Using the parameterization ( |(J.1|) for the ghost propagator, performing the color summa¬ 
tions and retaining only the part of T^g^ that is antisymmetric in the color indices one 
obtains from (|C.3|) 

' ^ t ^ dXsiniXK^) f dX'e-^^ (C.4) 

Jo Jo 


= ^^Q;(25 - i)d 


J (27r)^ 

where we have already taken the limit H —>■ 4 since the integral ( |C.4|) is UV and IR finite. 
The vertex in the limit of vanishing momentum transfer 0 is particularly simple to 

evaluate 




OdA 


(g^-O) = 


-abc-c 

-- 7 


a(26 — 1)6^ sm(K^X) 




dX- 


- 1)P 


A 


327rK^ 


(C.5) 


Using that a{26 — 1 ) = 3 in CCG, the last expression in ( |C.5|) for the antisymmetric 
vertex at vanishing momentum transfer becomes ( b.l 8 |) of the main text. 


D RG-invariance of the effective long-range interac¬ 
tion 


We here determine the renormalization of the effective interaction (iH) to order for 
planar contributions to large Wilson loops of the kind shown in Fig. 2. Power counting 
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suggests that UV-divergences do not arise from this effective interaction since the coupling 
K has positive mass dimension. To lowest order in g{^) the interaction (|6.7|) is just the 
product of the gluon propagator and the coupling g^. This interaction therefore depends 
on the renormalization point fi. We will show that this /i-dependence of the lowest order 
interaction is however canceled to order g"^ by vertex corrections in the Wilson loop. 


The Wilson loop itself is a RG-invariant [p0[| observable, and can therefore be expressed 
in terms of the bare coupling gs and the bare gluonic held Ab which do not depend on 
the renormalization scale. The bare held and coupling are related to the corresponding 
renormalized quantities by ( p.l|) , with the renormalization constants of (|3.2|) . The product 
of propagator and coupling constant g"^ renormalizes with Z'^Z^ 

gl {AB{x)ABm = Z^Z.g^ (/l(x)R(0)) (D.l) 

(p.l|) is proportional to the ehective 2-point interaction (|6.3|) in the Wilson loop however 
only to lowest order in g^. In order g'^ the divergent parts of the diagrams in Fig. 6 
also contribute to the ehective 2-point interaction (|6.7|) . As we will see, these divergences 
compensate the factor ZgZ^ in (|D.1[) to order g^. This implies that the (hnite) corrections 
from the contributions of Fig. 6 precisely cancel the /i dependence of the renormalized 
ehective 2-point interaction to order g"^. The anomalous dimension of the ehective 2- 
point interaction (|6.7|) is then at least of order g'^, which is sufficient to prove that it is 
renormalization group invariant in the limit fi ^ oo. 


We show that the divergent corrections from planar diagrams of order g^ shown in 
Fig. 6 can be absorbed in the longitudinal interaction by a multiplicative factor 

(ZlZ^W = l + ^d^^[3 + a) + 0{g*) (D.2) 

if the lowest order interaction is longitudinal (i.e. we consider only 0{g^) corrections to 
the long-range interaction). 


Only the self-energy insertions (Fig. 6a) and vertex correction (Fig. 6b) are UV- 
divergent in dimensional regularization as £ = {4 — D)/2 ^ 0. To evaluate the diver¬ 
gences, only the behavior of propagators at high momenta are needed and perturbative 
propagators and vertices are therefore sufficient in an asymptotically free model. 



Fig. 6: Planar corrections of order to the effective longitudinal interaction. Sections of 
the Wilson loop are shown in bold and the path ordering is indicated by arrows; the effective 
longitudinal interaction is depicted as a thin line and pertnrbative glnon exchanges as wavy 
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lines, a) A self-energy insertion on the Wilson line, h) K perturbative vertex corrections to the 
long-range longitudinal gluon interaction, c) Counter-term Z‘^Z^ — 1 due to renormalization of 
the effective longitudinal interaction, d) A. finite perturbative correction of order g^. 


Let us first analyze the vertex correction of Fig. 6b). In the Wilson loop this correction 
gives rise to a loop integral of the form 




Day{p + k) (D.3) 


(27r 

xD^p{k) \6p^{k‘^ -{p+ kf) - kpk^ + (p + k)p{p -b k). 


where we have used that the lowest order effective interaction is longitudinal. In 
the color traces have been performed, Dpu{p) denotes the perturbative gluon propagator 
and x{9),y{9') are points on the Wilson loop. The ©-function in ( p.3|) arises due to path 


ordering. Note that the term in square brackets of ( p.3|) is a linear superposition of trans¬ 
verse projectors. Inserting the explicit expression for the perturbative gluon propagator, 
the integrand of (S) can be rewritten as 


DfMp{k) [6p^{k^ - (p + kf) - kpk„ + (p + k)p{p -b k)a\ D„fp + k) = D^fp + k)- 

n , kpp^+pp{p + k)^ , ^ IN/, p,2 , fk'^ + pk){p‘^ + 2pk) 

—Dpv\k) H-— 1 \2i 2 -b (a — l)(fc^(p -b k)v{k -b pk)- 


(p -b k)‘^kP‘ 


k^{p + ky 


(D.4) 


The hrst two terms in the last expression for the integrand of (|D.3|) precisely compensate 
the divergent contributions from the self-energy insertions^® of Fig. 6a). The other terms 
in (p.4|) are longitudinal in one or the other (or both) gluon vertices. The /i-dependent 
divergent contributions of order from the vertex- and self-energy corrections to the 
effective interaction therefore arise from the divergent part of the integral 


PJ^I{P)\ 


Fig.Qa-\-b 

d^k 


convergent + [ j d9 dff Q{9 - 9’)^'^x 

2 J J d9 d9' 

1 ( d d \ 


X 


+ (o 


(27r)^ 

- 1 ) 


(A; — p/2)2(/c-b p/2)2 
pkijy — p^/4) 


-P^ 

d 


dx. 


{k + p/2)^{k — p/2)"^ 
(D.5) 


{k-p/2)p- -(A; + p/2), 

oyu 


_d_ 

dx,, 


^i[{k-p/2)x-(k+p/2)y] 


Observing that the remaining integral is divergent only for x = y, one can now perform 
one of the line integrals and obtains 


Pc.r“(p)| 


Fig.6a-\-b 


= convergent -b g^nt" 


d9 


PpX 


^®Note that the long-range interaction in Fig. 6a) is longitudinal - /r-dependent divergences arise only 
from so-called “pinch” contributions ^l| with at least two coinciding vertices. Apart from a difference in 
the ove rall sign, the integrands of these “pinched” self-energy contributions are just the first two terms 
in (Q. 
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X 


d^k 

(27r)^ 


{k — p/2y{k + p/2)2 


+ 


U ‘^^TL 

= convergent + , ^ —(3 + a)t“ 

QAir^e 


dO 


a — 1 k‘^{k‘^ — p^/4:) 

D {k + p/2Y{k-p/2Y 


dx 


de 




(D.6) 


as the divergent part of the vertex and self-energy corrections to the effective longitudinal 
interaction in order g^. The divergent contributions we calculated arise at both vertices 
of the effective interaction and thus precisely cancel the counter-term — 1 from the 
renormalization of the lowest order effective interaction to order g'^. In the Wilson loop, 
the p-dependence of the effective longitudinal interaction is therefore compensated by the 
//-dependence of vertex and self-energy corrections to order g"^. 


Only the perturbative diagrams considered here are planar and UV-divergent to order 
g^. It might seem that IR-divergent contributions arise from the exchange of longitudinal 
perturbative gluons shown in Fig. 6d). This is however not the case and the IR-divergences 
in fact cancel. To order g"^ the effective longitudinal 2-point interaction is independent of 
fi and its anomalous dimension is therefore of order g^. Perturbative corrections to the 
effective long range interaction are therefore analytic in g and K and s in (|6.7| ) approach 
RG-invariant constants as 0. 


E Integral equation for a planar Wilson loop with 
longitudinal interaction 


In the following we derive the integral equation ( |6.8| ) satished by a Wilson loop with an 
effective longitudinal gluonic 2-point interaction (|6.7|) . As has been argued in section 6, 
we need to consider only contributions which are of leading order in an 1/n expansion - 
for color-vector interactions, these are planar diagrams with all exchanges inside the loop. 

Consider a typical contribution to the Wilson loop with only 2-point interactions of 
this type. We open the loop at some point A as indicated in Fig. 7. Since the vertices 
on the loop are ordered and the exchanges are planar in the above sense, there is a one- 
to-one correspondence between one-sided planar contractions of vertices on a “line” and 
corresponding contributions to the Wilson loop if each topologically different one-sided 
planar diagram of the “line” is counted exactly once. Although several topologically 
different diagrams on the “line” may correspond to the same diagram for the Wilson 
loop, this multiplicity precisely accounts for the statistical factor of that diagram in the 
loop. We are thus left with the combinatorial problem of generating all topologically 
different one-sided planar contributions to the “line” exactly once. 

^°Note that for a transverse gluonic interaction the divergent perturbative corrections could not be 
absorbed multiplicatively, since the operators and [Api,Av\ mix in order g^. 
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A 

Fig. 7; A planar contribution to the Wilson loop from a 2-point interaction and the corresponding 
one-sided planar “line”-diagram when the loop is opened at A 

The closed curve of the Wilson loop C := : 6 G [0, 27r] with = x^(27r)}, 

is parameterized by 0 G [0, 27r]. In the following we reserve Greek letters for the one¬ 
dimensional coordinates on the loop respectively “line”. An effective longitudinal inter¬ 
action of the form — y) = 5°'’^d^dya{{x — yY)/n in leading order of 1/n results in a 
scalar factor 

P) = (E-1) 

for each interaction on the “line” between points a and (3 of such a one-sided planar 
diagram, since the color matrices at each vertex always contract to give the quadratic 
Casimir of the fundamental representation. In leading order of 1/n, the “line” is therefore 
proportional to the unit matrix. Note that for a transverse exchange the interaction 
depends on the direction relative to the line element and is not of the form ( [E.l| ). For 
an abelian group there is no distinction between the ordering of the vertices and any 
longitudinal interaction gives a vanishing contribution to the Wilson loop, since ( [E.l|) is a 
total derivative in the loop parameter. In non-abelian gauge theories, the ordering of the 
vertices does however matter and longitudinal interactions do contribute to the Wilson 
loop. Without the spontaneous symmetry breaking, the contribution from longitudinal 2- 
point interactions would be compensated by contributions from Green’s functions of higher 
order. As explained in the main text, we however do not expect such a compensation of 
the long range longitudinal interaction when the global SU (n)-symmetry is spontaneously 
broken. 

Let us denote by Fk{a,l3) the sum of all one-sided planar diagrams containing k 
longitudinal interactions between the points a and (3 on the “line”, with each topologically 
distinct diagram counted exactly once. Diagrammatically, Ffc+i(a, (3) is related to Fi{a, (3) 
with I < k by: 





Fig. 8; Iterative equation for Fk+i{a, f3). 

Fig. 9 expresses the fact that the hrst interaction in a one-sided planar diagram of order 
k + 1 “encloses” a one-sided planar diagram of order less than k and that the remaining 
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diagram after this first interaction is one-sided planar and of the complimentary order. 
Naturally, one has to sum over all intermediate orders and integrate over all possible 
(ordered) positions of the vertices for the hrst interaction between the points a and (5. It 
is not hard to show inductively that each topologically distinct one-sided planar diagram 
of order k + 1 appears exactly once in Ffc+i. Written explicitly, the one-sided planar 
contribution of order k + 1 therefore is, 

Ffc+i(a,/?)=^ f d-f f dru(r,7)F^(r,7)Ffc_^(7,/?) (E.2) 

r» ^O' Jo 

m=0 

where Fo{a, f3) = 1 is the contribution with no interactions. The total one-sided planar 
contributions W{a,(3) to the “line” between a and /? > a is just the sum of all F^, since 
each has the statistical factor 1: 


W(a, P) = Y1 FJa, (5) = Fo{a, /3) + ^ F^+M, (3) (E.3) 

fc =0 fc =0 

Using ( [E.2|) and rearranging the double sum, one obtains 

^ ^ rjd p(3 

W{a,f3) = Fo{a, (3) + / dTv{T,'y)Fm{T,'y)Fk-m{'y, (d) 

fc= 0 m= 0 ‘^“ 

= 1+ f dj f dr u(r, 7 )lU(r, 7 ) 1 U( 7 ,/?) 

Jo Jo 

= 1+ f dr f d 7 u(r, 7 )lU(r,7)1U(7,/9). (E.4) 

Jo Jr 


Relatively simple is the case of a circular Wilson loop parameterized by an angle with 
a longitudinal interaction of the form (pll.ll) . The distance between points on the circular 
loop in this case depends only on the difference in their angular positions and the radius 
R-. {x — uY = 4i?^sin^ Eor a longitudinal and translationally and rotationally 

invariant interaction, u(a,/d) = VR{a — f3) of o in this case is a function of the angular 
difference a — (3 and the radius only. Thus, for a circle of radius R, W{a, (3) = — (3) 

also depends only on the angular difference and one can rewrite (|E.4|) as: 


WR{e) 


1 + d9' f dF'vRie" - e')WR{e" - e')WRie - e") 

Jo Jd' 

1 + de' f de"vR{e")WR{e")WR{e'-6"). 

Jo Jo 


(E.5) 


To solve (|E.5| ) numerically, it is advantageous to consider instead the equivalent non-linear 
integro-differential equation: 


= r de'VR{e')WR{e')WR{e - e'), 

ad Jo 


(E.6) 


with the initial value 1 Ur( 0) = 1. The expectation value of the Wilson loop in this 
approximation is then nWR{27r). 
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